arXiv:1501.03679vl [math.PR] 15 Jan 2015 


KINETIC BROWNIAN MOTION ON RIEMANNIAN MANIFOLDS 


J. ANGST, I. BAILLEUL, AND C. TARDIF 


Abstract. We consider in this work a one parameter family of hypoelliptic diffusion pro¬ 
cesses on the unit tangent bundle of a Riemannian manifold (A4,g), collectively called 

kinetic Brownian motions, that are random perturbations of the geodesic flow, with a pa¬ 
rameter CT quantifying the size of the noise. Projection on A4 of these processes provides 
random paths in A4. We show, both qualitively and quantitatively, that the laws of these 
M-valued paths provide an interpolation between geodesic and Brownian motions. This qual¬ 
itative description of kinetic Brownian motion as the parameter a varies is complemented by 
a thourough study of its long time asymptotic behaviour on rotationally invariant manifolds, 
when cr is fixed, as we are able to give a complete description of its Poisson boundary in 
geometric terms. 
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1. Introduction 

1.1. Motivations and related works. We introduce in this work a one parameter family 
of diffusion processes that model physical phenomena with a hnite speed of propagation, 
collectively called kinetic Brownian motion. The need for such models in applied sciences 
is real, and ranges from molecular biology, to industrial laser applications, see for instance the 
works |CDR08l IGS13[ IGKSW14] and the references therein. As a hrst step in this direction, 
we consider here what may be one of the simplest example of such a process and provide a 
detailed study of its behaviour in a fairly general geometric setting. In the Euclidean space 
kinetic Brownian motion with parameter a, is simply described as a random path {xt)t>o 
with Brownian velocity on the unit sphere, run at speed so 

( 1 - 1 ) ^ = xt, xt = 

at 

for some Brownian motion W on the unit sphere of M'^. 



Figure 1. Kinetic Brownian motion in Euclidean space, velocity on the sphere 
on the left and the corresponding position on the right. 

In contrast with Langevin process, whose RAyalued part can go arbitrarily far in an arbi¬ 
trarily small amount of time, kinetic Brownian motion provides a bona hde model of random 
process with hnite speed. Its dehnition on a Riemannian manifold follows the intuition pro¬ 
vided by its version, and can be obtained by rolling on A4 without slipping its Euclidean 
counterpart. Figure 2 below illustrates the dynamics of kinetic Brownian motion on the torus, 
as time goes on. 





Figure 2. Simulations of a the kinetic Brownian motion on the torus over 
different time intervals. 

We devote most of our efforts in this work in relating the large noise and large time behaviour 
of the process to the geometry of the manifold. On the one hand, we show that kinetic 
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Brownian motion interpolates between geodesic and Brownian motions, as a ranges from 0 to 
oo, leading to a kind of homogenization. Onr nse of rongh paths theory for proving that fact 
may be of independent interest. We hrst prove the interpolation property in the model space 
and strengthen the associated convergence resnlts into some rongh paths convergence 
resnlts. The twist here is that once the latter result is proved, the fact that kinetic Brownian 
motion can be constructed from the rough path lift of kinetic Brownian motion in by 
solving a rough differential equation, together with the continuity of the Ito map, in rough 
path topology, provides a clean justification of the homogenization phenomenon. 

This result strongly echoes Bismut’s corresponding result for his hypoelliptic Laplacian 
[BisOh] which, roughly speaking, corresponds in its simplest features to replacing the Brownian 
velocity on the sphere by a velocity process given by an Ornstein-Uhlenbeck process. As a 
matter of fact, our method for proving the above mentionned homogenization result can also 
be used to recover the corresponding result for Bismut’s hypoelliptic diffusion. 


On the other hand, we are able to give a complete description of the Poisson boundary of 
kinetic Brownian motion when the underlying Riemannian manifold is sufficiently symmetric 
and a is fixed. This is far from obvious as kinetic Brownian motion is a hypoelliptic diffusion 
which is non-subelliptic. We take advantage in this task of the powerful devissage method 
that was introduced recently in |AT14j as a tool for the analysis of the Poisson boundaries of 
Markov processes on manifolds. Its typical range of application involves a diffusion {zt)t>o that 
admits a subdiffusion {xt)t>o whose Poisson boundary is known. If the remaining piece yt of 
Zt = {xt, yt) converges to some random variable yoo, the devissage method provides conditions 
that garantee that the invariant sigma held of {zt)t>o will be generated by y^o together with the 
invariant sigma held of {xt)t> 0 j see the end of Section 3.2.2 where the main results of |AT14] 
are recalled. In the present situation, and somewhat like Brownian motion on model spaces, 
the Poisson boundary of kinetic Brownian motion is described by the asymptotic direction in 
which the process goes to inhnity. It is remarkable, however, that depending on the geometry, 
kinetic Brownian motion may have a trivial Poisson boundary while Brownian motion will 
have a non-trivial Poisson boundary. 


Kinetic Brownian motion is the Riemannian analogue of a class of dihusion processes on 
Lorentzian manifolds that was introduced by Franchi and Le Jan in |FLJn7] . as a generalization 
to a curved setting of a process introduced by Dudley |Dnd66] in Minkowski spacetime. These 
processes model the motion in spacetime of a massive object subject to Brownian huctuations 
of its velocity. Despite their formal similarities, the causal structure of spacetime makes 
relativistic dihusions very diherent from kinetic Brownian motion, as the reader will hnd out 
by reading the litterature on the subject, such as |FLJ07l IFLJlOt IFraOQt IBaiOSt IBRlOt IBailOt 
Angl4|, for instance. 


We have organized the article as follows. Kineti c Bro wnian motion on a given Riemannian 
manifold {Ai,g) is introduced formally in Section [L2 below. Like its elementary RAyalued 
version, it lives in the unit tangent bundle of AJ, and it has almost surely an inhnite lifetime 
if the manifold is geodesically complete. Section is dedicated to proving that the manifold¬ 
valued part of a time rescaled version of kinetic Brownian motion converges weakly to Brownian 
motion as the intensity of the noise, quantihed by increases indehnitely, if the manifold 
is stochastically complete. This is the main content of Theorem 2.1.2[ which is proved using 
rough paths theory. The necessary material on this subject is recalled, so the reader can follow 
the proof without preliminary knowledge about rough paths. Section provides a thourough 
description of the asymptotic behaviour of kinetic Brownian motion on a rotationnaly invariant 
manifold, through the identification of its Poisson boundary in a generic setting. 
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We collect here a number of notations that will be used throughout that work. 


• We shall use Einstein’s well-known summation convention and (A4,g) will denote a 

d-dimensional oriented complete Riemannian manifold, whose unit tangent bundle and 
orthonormal frame bundle will be denoted respectively by and OAi. We shall 

denote hy z = (x, e) a generic point of OAi, with x E At and e : —)■ T^Ai, an 

orthonormal frame of T^Ai] we write tt : OA4 -E At for the canonical projection 
map. Last, we shall denote by (ei,..., e^) the canonical basis of M'’*, with dual basis 

(el, ■■■,4)- 

• Denote by (Vi) 2 <i<d, the canonical vertical vector helds on OAl, associated with the 

Lie elements Vi := 0 — €i 0 e* of the orthonormal group of The Levi-Civita 

connection on TAl dehnes a unique horizontal vector held ifi on OAi such that 
(7r*iLi)(2;) = e(ei), for all 2 ; = (x, e) E OAi. The how of this vector held is the natural 
lift to OAi of the geodesic flow. Taking local coordinates x* on Ai induces canonical 
coordinates on OAi by writing 

d 

Ci := e{ei) = ^eld^j. 
i=i 


Denoting by TL the Christohel symbol of the Levi-Civita connection associated with 
the above coordinates, the vector helds Vi and Hi have the following expressions in 
these local coordinates 


VM 


Hdz) 


= ef 


= e 


A 

d 

^dxi 


^ de'l ’ 
! 

ij\ 


j ^ 
de’l ■ 


2 <i < d, 


1.2. Definition of kinetic Brownian motion. As said above, kinetic Brownian motion on 
a d-dimensional oriented complete Riemannian manifold (Ai,g), is a dihusion with values in 
the unit tangent bundle T^Ai of Ai. In the model setting of M'’*, it takes values in x 
and is described as a random path run at unit speed, with Brownian velocity, as described 
in Equation (1.1). As in the classical construction of Eells-Elworthy-Malliavin of Brownian 
motion on At, it will be convenient later to describe the dynamics of kinetic Brownian motion 
on a general Riemannian manifold {Ai,g) as the projection in T^At of a dynamics with values 
in the orthonormal frame bundle OAi of At, obtained by rolling without splitting kinetic 
Brownian motion in The following direct dynamical dehnition in terms of stochastic 
diherential equation provides another description of the dynamics of kinetic Brownian motion 
which we adopt as a dehnition. The equivalence of the two point of views is shown in section 
2.4. 1[ In the sequel a will always stand for some non-negative constant which will quantify 


the strength of the noise in the dynamics of kinetic Brownian motion. 


Definition 1.2.1. Given zq E OAi, the kinetic Brownian motion with parameter a, started 
from zq, is the solution to the OAi-valued stochastic differential equation in Stratonovich form 

(1.2) dzt = Hflzt) dt + aVflzt) odBl, 

started from zq. It is defined a priori up to its explosion time and has generator 

i=2 
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It is elementary to see that its canonical projection on is a diffusion on its own, also 

called kinetic Brownian motion. In the coordinate system := (x\e\) on T^Ai induced 

by a local coordinate system on Ai, kinetic Brownian motion satishes the following stochastic 
differential equation in ltd form 


(1.3) 


dxl = x\dt, 

2 

dx\ = dt + adMl —^{d — l)x] dt. 


where 1 < i < d, and where Mt is an M'^-valued local martingale with bracket 


d{M\M^)t = (^g'^{xt) -xlxl^dt, 


for any 1 < i, j < d and stand for the inverse of the matrix of the metric in the coordinates 
used here. 


The readers acquainted with the litterature on relativistic diffusions will recognize in equa¬ 
tion (1.2) the direct Riemannian analogue of the stochastic differential equation dehning rel¬ 
ativistic diffusions in a Lorentzian setting. Despite this formal similarity, the two families of 
processes have very different behaviours. As a trivial hint that the two situations may differ 
radically, note that the unit (upper half) sphere EI'^ in the model space of Minkowski 
spacetime, is unbounded. As a result, there exists deterministic HAvalued paths Xg that ex¬ 
plode in a hnite time, giving birth to exploding M^’Avalued paths Xg = xq + J^Xrdr. The 
work jBai m even gives some reasonnable geometric conditions ensuring the non-stochastic 
completeness of relativistic diffusions. No such phenomenon can happen in or on a com¬ 
plete Riemannian manifold for a path dehned on a hnite open interval, and run at unit 
speed. 


Proposition 1.2.2 (Non-explosion). Assume the Riemannian manifold {Ai,g) is complete. 
Then kinetic Brownian motion has almost surely an infinite lifetime. 


Proof. Denote by r the lifetime of kinetic Brownian motion zt = {xt,et), and assume, by 
contradiction, that r is hnite with positive probability. Since the path {xt)o<t<T has unit 
speed, it would converge as t tends to r, on the event {r < cx)}, as a consequence of the 
completeness assumption on Ai. The horizontal lift ■= {xt,et)o<t<T of {xt)o<t<T in 

OAi would converge as well. Write Vi for the Lie element ® ej — ei ® e* of the orthonormal 
group SO{d) of M'^. We have then e* =: Ctht where the process (h*) G SO{d) satishes the 
stochastic diherential equation 

dht = —Viht odBfi 

in particular (hj is well dehned for all times t > 0. Consequently, both processes (cj) and (zt) 
should converge as t tends to r, contradicting the necessary explosion oi Zt . □ 


From now on we shall assume that the Riemannian manifold {A4,g) is complete, and turn 
in the next section to the study of kinetic Brownian motion as a function of a. 


2. From geodesics to Brownian paths 

Let emphasize the dependence of kinetic Brownian motion on the parameter a by denoting 
it zf = {xf,ef) e OAi. We show in this section that the family of laws of x^ provides a kind 
of interpolation between geodesic and Brownian motions, as a ranges from zero to inhnity, 
as expressed in Theorem 2.1.2 below and illustrated in Fig. below when the underlying 
manifold is the 2-dimensional hat torus. 
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Figure 3. Kinetic Brownian motion on the flat torus for cr = 10 10 1,2,4,10. 

2.1. Statement of the results. To fix the setting, add a cemetary point d io M., and endow 
the union with its usual one-point compactiflcation topology. That being done, denote 

by flo the set of continuous paths 7 : [0, 1] —)• Af U {5}, that start at some reference point xq 
and that stay at point d if they exit the manifold M.. Let T := Vte[o 1] where (JT)o<t<i 
stands for the filtration generated by the canonical coordinate process. Denote by Sr the 
geodesic open ball with center x<^ and radius i?, for any R> 0. The first exit time from Sr is 
denoted by tr, and used to define a measurable map 

Tr : flo ^ ^*([0,1], Sr), 

which associates to any path ( 74 )o<t<i G Dq the path which coincides with 7 on the time 
interval [O, tr] , and which is constant, equal to 7 ^-^, on the time interval [tr, l]. The following 
definition then provides a convenient setting for dealing with sequences of random process 
whose limit may explode. 

Definition 2 . 1 . 1 . A sequence (Pn)n>o of probability measures on is said to converge 

locally weakly to some limit probability P on (p,o,iF) if the sequence P„ o Tr^ of probability 
measures on C([0,1], Sr) converges weakly to F o , for every S > 0. 

Equipped with this definition, we can give a precise sense to the above interpolation between 
geodesic and Brownian motions provided by kinetic Brownian motion. 

Theorem 2 . 1.2 (Interpolation). Assume the Riemannian manifold (A4,g) is complete. Given 
A = ( 2^05 Co) € OM. we have the two following asymptotics behaviours. 

• The law of the rescaled process converges locally weakly under F^^ to Brow¬ 
nian motion on M., run at speed time interval [0,1], as a goes to infinity. 

• The law of the non-rescaled process converges locally weakly under F^q, as 

a goes to zero, to a Dirac mass on the geodesic started from Xq in the direction of the 
first vector of the basis cq. 
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Implicit in the above statement concerning the small noise asymptotics is the fact that the 
geodesic curve is only run over the time interval [ 0 , 1 ]; no stochastic completeness assumption 
is made above. The proof of the second item in the interpolation theorem is trivial, since 
the generator of the process then converges to the generator of the geodesic flow as it is clear 
from the Definition 1.2.1 of the process. We shall hrst prove the result stated in the hrst 
item by elementary means in the model Euclidean case, in Section |2.2[ Using the tools of 


rough paths analysis, we shall see in Section 2.3 that elementary moment estimates allow 


a strengthening of the weak convergence result of Section 2^ into the weak convergence of 
the rough path lift of the Euclidean kinetic Brownian motion to the (Stratonovich) Brownian 
rough path. To link kinetic Brownian motion in Euclidean space to its Riemannian analogue, 
we use the fact that the latter can be constructed by rolling on M without slipping the former, 
using Cartan’s development map. This means, from a practical point of view, that one can 
construct the Al-valued part of kinetic Brownian motion as the solution of an Wl-valued 
controlled differential equation equation in which the Euclidean kinetic Brownian motion plays 
the role of the control. This key fact will enable us to use the continuity of the Ito-Lyons map 


associated with Cartan’s development map, and transfer in Section [2i4| the weak convergence 
result proved for the Euclidean kinetic Brownian motion to the curved setting of any complete 
Riemannian manifold. 

Note that the above theorem only involves local weak convergence; it can be strengthened 
under a very mild and essentially optimal natural assumption. 

Corollary 2.1.3. If the manifold (M,g) is complete and stochastically complete then the 
rescaled process (a^CT 2 t)o<t<i coniierg'es in law under as a goes to infinity, to Brownian 
motion run at speed over the time interval [0,1]. 


X.-M. Li proved recently in |Lil5] an interpolation theorem similar to Corollary 2.1.3 under 
stronger geometric assumptions on the base manifold M., requiring a positive injetivity radius 
and a control on the norm of the Hessian of the distance function on some geodesic ball. Her 
proof rests on a formulation of the weak convergence result in terms of a martingale problem, 
builds on ideas from homogenization theory, and uses tightness techniques. It is likely that the 
very robust nature of our proof, based on the rough paths machinery, offers a convenient setting 
for proving more general homogenization results at a low cost. As a basic illustration, notice 
for example that the proof below works verbatim with the Levi-Civita connection replaced by 
any other affine metric preserving connection H : TAl —?■ TOM.. The limit process in OM is 
not in that case the lift to OM of a Brownian motion on M anymore, but it is still described 
as the solution to the Stratonovich differential equation 

det = H(et) odBt, 

where B is an pAvalued Brownian motion. 

Note that a slightly more general family of diffusions on the frame bundle than those given 


by equation (1.2) is considered in |Lil5j . where in addition to the Brownian noise in the 
canonical vertical directions, scaled by a factor cr, she also considers a vertical constant drift 
independent of a. It is elementary to adapt our method to this setting. 

Remark 2.1.4. The idea of using rough paths theory for proving elementary homogenization 


results as in theorem 2.1.2 was first tested in the work |E(1L14] of Friz, Gassiat and Lyons, 
in their study of the so-called physical Brownian motion in a magnetic held. That random 
process is described as a path {xt)o<t<i in modeling the motion of an object of mass m. 
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with momentum p. = mx., subject to a damping force and a magnetic field. Its momentum 
satisfies a stochastic differential equation of Ornstein- Uhlenbeck form 

dpt = Mptdt + dBt, 
m 

for some matrix M whose eigenvalues all have positive real parts, and B is a d-dimensional 
Brownian motion. While the process {Mxt)Q<t<i is easily seen to converge to a Brownian 
motion, its rough path lift is shown to converge in a rough paths sense in h'^, for any q > 2, 
to a random rough path different from the Brownian rough path. 


The proof of the interpolation theorem 2.1.2 and Corollary 2.1.3 is split into three steps, 


performed in t he next three subsections. We prove the Euclidean version of Corollary 2.1.3 
in Section 2.2, and strengthen that weak convergence result in (C ([0,1],]R‘^) , || ■ ||oo) into a 
weak convergence result in rough paths topology of the rough paths lift of Euclidean kinetic 
Brownian motion. This is done in Section [2^ by using some general compactness criterion on 
distributions in rough paths space. The point here is that kinetic Brownian motion on any 
complete Riemannian manifold can be constructed from its Euclidean analogue by solving a 
controlled differential equation in which the control is the Euclidean kinetic Brownian motion. 
One can then use Lyons’ universal limit theorem, on the continuity of the Ito-Lyons map, 
to transfer the weak convergence result of the rough kinetic Brownian motion in M'’* to its 


Riemannian analogue; this is explained in Section 2.4, where Corollary 2.1.3 is also proved. 


2.2. Proof of the interpolation result in the Euclidean setting. When the underlying 
manifold M. is the Euclidean space M'^, the state space of kinetic Brownian motion becomes 

X inherited from the canonical 


^ gd 1 ^ with coordinates {x,x) = {x\x^)i<i<d € 
coordinates on M'’* x M'’*. System (1.3) describing the dynamics of kinetic Brownian motion in 


a general local coordinate system takes in the present setting the simple form 

dxl = x\dt, 


( 2 . 1 ) 


dxl = 


,d — 1 


-a 


• • 7 


)ml, 


-x‘,dt + {S‘^ 

for 1 < i < d, with W = (W ^,..., IT'’*) a standard M'’*-valued Brownian motion. 


Proposition 2.2.1. Given (a;o,xo) ^ \ we have the following two asymptotic regimes, 

in terms of a. 

(1) The non-rescaled process (a;^)o<t<i converges weakly to the Dirac mass on the geodesic 
path [xq + txo)^^_^^.^, as a tends to zero. 

(2) The time-rescaled process {x^ 2 f)o<t<i converges weakly to a Euclidean Brownian motion 

with covariance matrix times the identity, as a tends to infinity. 


Proof. The convergence result when a tends to zero is straightforward. We present two proofs 
of the weak convergence to Brownian motion as a tends to infinity, to highlight the elementary 
nature of this claim. In order to simplify the expressions, let us define for all t > 0 


Xf := 


The first approach takes as a starting point the integrated version of equation (2.1), namely 


( 2 . 2 ) 


Xf — Xo -f 















KINETIC BROWNIAN MOTION ON RIEMANNIAN MANIFOLDS 


9 


where is a d—dimensional martingale whose bracket is given by 




{M”\ M->), = i I (0 - ifi0) ds. 


The mid-terms in the right hand side of Equation (2.2) clearly go to zero when a tends to oo, 
uniformly in t, so that the asymptotic behavior of is the same as that of the martingale 
. To analyse that martingale, note that the time-rescaled process {yt)t>o ■= (it/o-2)t>o is a 
standard Brownian motion on solution of the equation, for each 1 < i < d 


dyl = 


d-i 


(2.3) 


2 

d-1 


yldt + Y .(0 -ylyi)dBi 

i=i 


-yldt -F dNl 


for some M'^-valued Brownian motion B. The bracket of the martingale is simply given in 
terms of {yt)t>o by the formula 




= {0-y^i)ds, 

so, for a hxed time t, the ergodic theorem satished by the process {yt)t>o entails the almost 
sure convergence 


lim {M", M-d), = ^ 

(T^ + OO ClyCl Ij 




The result of Proposition 2.2.1 then follows from the asymptotic version of Knight Theorem, 
as stated form instance in Theorem 2.3 and Corollary 2.4, pp. 524-525, in the book |RY99j of 
Revuz and Yor. The second approach consists in starting from the integral representation 




(2.4) 


— ^0 H-F 


cr^ 


y\ds, 


where {yt)t>o is the above standard spherical Brownian motion, so that the result can alterna¬ 
tively be seen a consequence of a standard central limit theorem for ergodic diffusions applied 
to the ergodic process {yt)t>o, see for instance the reference |CCG12j . □ 


We shall prove in Section 2.3 that the weak convergence result of Proposition 2.2.1 can be 


enhanced to the weak convergence of the rough path associated with the random path 
to the Stratonovich Brownian rough path. This requires the study of the 
valued process defined for any f > 0 by the integral 






For those readers not acquainted with tensor products, one can simply see the tensor product 


a (g) & of two vectors a,b in M“, as the linear map x G 


HA- 




Proposition 2.2.2. The -valued process (X^)i>o := converges weakly, 

as a goes to infinity, to the Brownian rough path over run at speed 
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Proof. It will be convenient to use the representation of Xf given by identity (2.4), for which 
there is no loss of generality in assuming xo = 0. Recall {Nt)t>o stands for the martingale part 

Q M'^-valr 

yf^ds + 


of ivt) i>o in Equation ( |2.3[ ), seen as an valued path. We have 

pcr'^t 


— _ 

* {d-l)a^ 


+ 


d-1^ 




(^(tH 


{d — 1 ) 2 ( 7 ^ 


—N^4^ 0 + 


dNs ®ys- 


dNs ® Ns 


As the symmetric part | of converges to the corresponding symmetric part of 

the Brownian rough path, by Proposition 2.2.1 , we are left with proving the corresponding 
convergence result for the anti-symmetric part of the process (X^)t>o, namely 




(2.5) 


(d-l)V4 
4 


yaH ® N^h - N^h < 8 ) yaH + 


dNs <S)ys-ys<^ dNs 


'-uH 


+ 


{d 


{Ns® dNs-dNs® Ns). 


As the process {yt)t>o lives on the unit sphere, the terms in the hrst line of expression (2.5) 
clearly converge to zero in L^, as a goes to infinity, uniformly on bounded intervals of time. 
Using the notations introduced in Equation (2.2), we are thus left with the martingale term 




H 


w / ® ®Ns)= / (m; 0 dMf - dMf 0 m;) . 

(^a — ij a Jq Jq 


The weak convergence of this term to the awaited Levy area of the Brownian rough path is 
dealt with by the following lemma, which concludes the proof. □ 


Lemma 2.2.3. The x -valued martingale (M^)t>o := (Mf, converges weakly 

as a goes to infinity, to the process 




{Bs 0 dBs 


dBs ®Bs) I 

t>o 


where B is an Mfi-valued Brownian motion run at speed 


Proof. Given x G M'^, let stand for the linear map from to x defined by the 

formula 


LxV := {v, vx* — XV *), 


where the star notation is used to denote linear form canonically associated with an element 
of the Euclidean space With this notation, we can write 

[ LM^dMf. 

Jo 

Recall that we have already proved that the process converges weakly when a goes to 
inhnity, to a Brownian motion Bt of variance . It follows from the continuity of the map 

X I— )■ Lx, that the process {Lm^, converges weakly to {Lb,, Refering to 

Theorem 6.2, p. 383 of |JS03j . the result of the lemma will follow from the previous fact if we 
can prove that the martingales are uniformly tight - see Dehnition 6.1, p. 277 of |JS03j . 
But since we are working with continuous martingales whose brackets converge almost surely 
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when a goes to infinity, we can nse Proposition 6.13 of |,TSn3] . p. 379, to obtain the awaited 
tightness. □ 


2.3. Prom paths to rough paths in the Euclidean setting. Proposition 2.2.2 shows that 
the natnral lift of as a weak geometric Holder p-rongh path, for 2 < p < 3, converges 
weakly to the Brownian rongh path, when seen as an element of a space of continnons paths 
on [0,1], endowed with the topology of uni form convergence. We show in this section that a 
stronger convergence resnlt holds, with the rongh path metric in place of the nniform topology. 


This stronger convergence resnlt will be the main ingredient nsed in Section 2.4 to prove the 
interpolation theorem [2.1.2 by relying on Lyons’ nniversal limit theorem. 

For the basics of rongh paths theory, we refer the reader to the very nice acconnt given in 
the book |FH14] by Friz and Hairer. Alternative pedagogical acconnts of the theory can be 
fonnd in |Bail4a[ IBanl3[ Lej09| , see also the book |FV10j of Friz and Victoir for a thonrongh 
acconnt of their approach of the theory. Let 2 < p < 3 be given. 

Proposition 2.3.1. The weak geometric Holder p-rough path (X^)o<t<i converges weakly as 
a rough path to the Brownian rough path. 

The remainder of this section is dedicated to the proof of this statement. Onr strategy of 
proof is simple. Using elementary moment estimates, we show that the family of laws of X°^ 
is tight in some rongh paths space. As the rongh path topology is stronger than the topology 


of nniform convergence on bounded intervals. Proposition |2.2.2| identifies the unique possible 
limit for these probability measures on the rough paths space, giving the convergence result 
as a consequence. 

To show tightness, we shall rely on the following Kolmogorov-Lamperti-type compactness 
criterion; see Corollary A.12 of |FV10j for a reference. As a shortcut, we will write Zts or Zt^s 
for the increment Zt — Zg, for any s < t, and any path with values in a vector space. 

Theorem 2.3.2 (Kolmogorov-Lamperti tightness criterion). Given any \ < 'y < consider 
the laws o/(X^)o<t<i = (X)’’, for a > 0, as probability measures on the metric space 

RP( 7 ) of weak geometric 'y-Holder rough paths. If the following moment estimates 


(2.6) 


sup E 

a 

sup E 


X 


ts 


X: 


ts 


^<Cq\t 
<Cg\t-S 


I £ 

s|U 


hold for all 0 < s < t < 1, for some positive constants Cq, for all q > 2, then the family of 
laws o/(X^)o<t<i is tight in RP( 7 ). 

We shall prove these two bounds separately, starting with the bound on Xf^. We use for 
that purpose the representation 

1 r"* 

= ^ / du, 

in terms of a unit speed Brownian motion {yt)t>o on the sphere, already used in the previous 
section, together with Equation (2.3) giving the dynamics of {yt)t>o in terms of an M'’*-valued 


Brownian motion Denoting by x* the coordinate of a vector x G M'’*, this gives the 

equation 

■ ' * AdRJ 


2 
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So we have 


(2.7) E 


X, 


ts 




I Va'^tjCr'^s \ 


E 


fcr^t 


- tM) dBi 


=:c,{(l) + (2)} 


for some positive constant Cq depending only on q. Note that term (2) contains an implicit snm 
over i and j. We use Burkholder-Davis-Gundy inequality to deal with it, taking into account 
the fact that the process {yt)t>o takes values in the unit sphere. This gives the estimate 


(2.8) 


(2) < (T“’’E 




- vivi) du 


Q , 1 9 

<22|t-s|T 


For term (1), remark simply that since the identity 


Vba 


d — 1 


yudu + Bba- / yu{yldBl) 


holds for all 0 < a < 6 < 1, for some valued Brownian motion S, we have 


< Cq 


|6 — a |'^ + |6 — a |2 + |6 — a |2 


using the fact that {yt)t>o lives on the unit sphere, so the process y^dW^ is a real-valued 
Brownian motion, together with Burkholder-Davis-Gundy inequality. So we have 

<cjb- a |2 

for all 0 < a < 6 < 1, for some well-chosen constant c > 2”?, since \yha\ < 2. The upper bound 
E[|Xts|'^] < Cq\t — s| 2 , follows then from equations (2.7) and ( |2.8 ). To deal with the double 
integral let us start from Equation (2.3) to write 

1 




(2.9) 






yvdv ®yudu 






I (T^g d 1 


yu^a^s G yudu 


{d - l)a^ y^4g 


0 y^du. 


Again, as the process {yt)t>o lives in the unit sphere, the first term in the right hand side above 
is bounded above by ^ \t — s|. Let us denote by (Ft)t>o the path in with values in 
whose velocity is given {yt)t>o- To deal with the term involving the martingale N in equation 
(2.9), we use an integration by parts, and the elementary inequality \a -|- < c (|a|'^ -|- |&|'^), 

to get the existence of a positive constant Cq depending only on q such that we have 

91 


( 2 . 10 ) 

where we have set 
(3) := E 


E 




A I ® VudU 


(7 


< Cg{(3) + (4)}. 




. 0-4*,(T^g 




91 


(4) := E 


1 

a"* 




— / dNu®Y^, 

a^s 


.du 


There is an absolute positive constant c for which 
(3) < cE 



9 

W.4t,cr4g 

1 

< cE 


2 q 

2 

E 


2 q 

cr^ 







cr^ 
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On the one hand, taking once more into account the fact the the process {yt)t>o lives on the 
unit sphere, the Burkholder-Davis-Gundy inequality gives us the upper bound 


E 






2q 


= E 


{S.j - y^^uVU^dBi 


2q 


< cAt-s\ 


On the other hand, we can use the identity 


(J^ 




and the fact that cr~‘^N^4^^^4g has the same law as — y^4^y^^4^_^)dBl, to see that 


E 


Y.4t, 




2q 


< < <! E 


ya'^tjCr'^s 




2q 


+ E 


( 5 '^ - ya*u'!/^4^)dBi 


2q 


We recognize in the first term on the right hand side term (1), with 2q in the role of g, and use 
the Burkholder-Davis-Gundy inequality to deal with the other expectation. The two results 
together give an upper bound of size |f — s|'^, up to a multiplicative constant depending only on 
q, showing that (3) is also bounded above by a constant multiple of |f — s|^. The inequality of 
Burkholder-Davis-Gundy is used once more, together with the bound just proved, to estimate 
term (4) by a constant multiple of 



Note that all the above estimates hold regardless of the values of a > 0. All together, these 
estimates prove that the two conditions of the compactness in Theorem 2.3.2 hold for the 
family of weak geometric y-Holder rough paths (X^)o<t<i. The family of laws of (X^)o<t<i is 
thus relatively compact on the space RP( 7 ). As its unique possible cluster point is identihed 
by Proposition |2. 2. 2] this proves the weak convergence of the random weak geometric 7 -Holder 


rough paths (X^)o<t<i to the Brownian rough path, as stated in Proposition 2.3.1 


2.4. Prom the Euclidean to the Riemannian setting via Cartan’s development map. 

The homogenization result proved in Proposition |2.3.1 puts us in a position to use the ma¬ 
chinery of rough differential equations and prove homogenization results for solutions of rough 
differential equations driven by X^^, using to our advantage the continuity of the ltd map in a 
rough paths setting. This is in particular the case of kinetic Brownian motion on any complete 
Riemannian manifold, which can be constructed from kinetic Brownian motion on using 
Gartan’s development map. The interpolation theorem 2 . 1.2 will follow from this picture of 
kinetic Brownian motion as a continuous image of X^^. Before following that plan, we recall 
the reader the basics of Gartan’s development method and rough differential equations. 
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2.4.1. Cartan’s development map. As advertized above, one can actually construct kinetic 
Brownian motion on a complete Riemannian manifold Ai by rolling on M. without slipping 
its Euclidean analogue. While this will be clear to the specialists from the very definition 
of kinetic Brownian motion, Cartan’s development procedure can be explained to the others 
as follows. In its classical form, this machinery provides a flexible and convenient way of 
describing paths on Ai from MAyalued paths; it requires the use of the frame bundle 
and the horizontal vector helds associated with a choice of connection (Levi-Civita connection 
presently). 

Let H stand for the TO Ad-valued 1-form on uniquely characterized by the property 

7r*(H(M))(2;) = e(n), 
for any u G and = (x, e) G OAi. 

Definition 2.4.1. Given zq = {xo,eo) in OAi, Cartan’s development of an M'’*-valued 
path (mt)o<t<i of class is defined as the solution to the ordinary differential equation 

( 2 . 11 ) dzt = il{zt) dmt 

started from zq. As before, we shall write Zt = {xt,et) G OAi. 


This description of an OAd-valued path may seem somewhat different from the kind of 
dynamics described by Equation (1.2) dehning kinetic Brownian motion. To make the link 


clear, assume, without loss of generality, that (mt)o<t<i is run at unit speed, and denote its 
speed by {rht)o<t<i. Reca ll the definition of the elements Vi of the Lie algebra of SO{d) given 
at the end of Section Then, given an orthonormal basis /o of M'’*, with /o(ei) = Aiq, solve 
the S'0((i)-valued ordinary differential equation 

dft = {ft{ei),dmt)vi{ft), 

started from /o, and define the M'^“^-valued control (ht)o<t<i, started from zero, by the formula 

dhl = {ft{ei),dmt), 

ioT 2 < i < d. Now, given any orthonormal basis go of Tm^Ai with go{ei) = eo(ei), it is 
elementary to see that the OAd-valued path obtained by parallel transport of go 

along the path (a;t)o<t<i = {'^{zt))o<t<i, satishes the ordinary differential equation 

(2.12) dz[ = Hi{z[)dt + Vi{z[)dhi. 

Let emphasize the fact that, by construction 

7l{zl) = Xt, 


Both Equations (2.11) and (2.12) make perfect sense with a path {mt)o<t<i whose derivative 


(mt)o<t<i is a continuous semimartingale; we talk in that case of stochastic development. 
The reader can consult the book |Hsu00j of Hsu for a pedagogical account of stochastic dif¬ 
ferential geometry. Considering in particular the stochastic development in Ai of kinetic 
Brownian motion on provides a control 

(2.13) dhl = a odBl, 


for some Stratonovich (d — l)-dimensional Brownian increment, in which case Equation ( |2.12[ ) 
is nothing but the equation giving the dynamics of kinetic Brownian motion, identifying 
with {z'f)t>o- Since the control given by formula (2.13) has a law independant of the arbitrary 
choice of frame fo satisfying the above constraint, the path x" = Tr(z.) = 7i{z'i') is seen to be 
the projection on Ad of a kinetic Brownian motion on Ai. 
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Although we know that the RAvalued part (m^)i>o of kinetic Brownian motion on con¬ 
verges weakly to some Brownian motion B run at speed , the setting of ltd or Stratonovich 
differential calculus is not robust enough to infer from that result the weak convergence of 
the developped process (a;^)o<t<i to the process (vr(wt))o<t<i, with {wt)o<t<i solution of the 
OAd-valued stochastic differential equation 

2 

dwt = H{wt) odBt, 

y a(a — 1) 

for which TT{wt) is nothing else than a scalar time-changed Brownian motion on M.. This 
comes from the fundamental lack of continuity of the solution map for stochastic differential 
equations, in ltd stochastic integration theory. This missing crucial continuity property is 
precisely what rough paths theory provides. We give in the next subsection the information 
on rough paths theory and rough differential equations needed to understand our reasoning, 
and refer the reader to the litterature on the subject for more insights; see the references 
below. 


2.4.2. Rough differential equations and the interpolation result. Let A/" be a smooth finite 
dimensional manifold. We adopt in this work the definition of a solution path to an A/'-valued 
rough differential equation given in jBail5| IBail4b] , as it is perfectly suited for our needs. It 
essentially amounts to requiring from a solution path that it satisfies some uniform Taylor- 
Euler expansion formulas, in the line of Davie’ seminal work |Davn7] . Let Ai,..., be vector 
fields of class on W; write A for the vector field-valued 1-form 

i 

A : M G HA- M*Aj. 

i=\ 

Let 2 < p < 3 and a weak geometric Holder p-rough path X = {X, X) over be given. 


Definition 2.4.2. An J\f -valued continuous path {mt)o<t<T is said to solve the rough differen¬ 
tial equation 

(2.14) dmt = A{mt)X{df) 

if there exists a constant a > 1 with the following property. For any 0 < s < t, there exists 
an open neighbourhood of rUg such that the Taylor-Euler expansion 

fim) = f{mg)+Xl{Aif){mg)+Xl^{AjAkf){mg) + 0{\t-sf) 

holds for all t close enough to s, for rut to belong to V*, and for any function f of class 
defined on V^. 


The remainder term 0(\t — s|“) is allowed to depend on X and /. It should be clear on this 
definition that the notions of classical controlled and rough differential equations coincide if 
the driving rough path is the canonical lift of a path, for the X^* term, of size {t — sfi in 
that case, can be put in the remainder. The crucial point for us here will be that if X is the 
(Stratonovich) Brownian rough path, associated with Brownian motion B, then this notion of 
solution gives back the solution of the Stratonovich differential equation 


dmt = Viimt) odBf 

See for instance the lecture notes |FH14j of Friz and Hairer, or the book LQ02 of Lyons and 
Qian. Classical results show that such a rough differential equation has a unique maximal 
solution started from any given point, and that the Ito map, that associates to the driving 
signal X the solution path x to the rough differential equation (2.14) is continuous in the 
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following sense. Fix T < t and cover the compact support of the path {mt)o<t<T by hnitely 
many local chart domains and (Oi)i<i<Ar, with O' C 0* for all 1 < i < iV. Then, 

there exists a positive constant 6 such that for every rough path Y 5-close to X in the Holder 
rough path distance, the solution path {yt)o<t<T to the rough differential equation 

dyt = A{yt)Y{dt) 

started from y^ = mQ, is well-dehned on the time interval [0,T] and remains in the open 
neighbourhood IJ^i of the support of {mt)o<t<T, with yt in Oi whenever m* is in O'. This 
very strong continuity property implies in particular the following “transport property” of 
weak convergence. We endow the space Hq dehned in the introduction of Section with the 
topology generated by the elementary open sets 

|( 2 /t)o<t<i e O([0,1],A/'U {9}) ; dist(i/t, 7 t) < e, for all t < r^jj, 

where ( 7 t)o<t<i is any element of Hq started from a point of the manifold, R and e are arbi¬ 
trary positive constants, and dist stands for any Riemannian distance function dehned in a 
neighbourhood of the support of 7 . The constant path, equal to d, is isolated. 


Proposition 2.4.3 (“Transporting” weak convergence by the ltd map). Let (X^)fc>o be a 
sequence of random weak geometric ^-Holder rough paths overM.^, whose distribution converges 
weakly in the rough paths space as k goes to infinity, to the distribution of some random weak 
geometric Holder p-rough path X. Fix xq G Af, and consider the solution path (m^)o<t<i to 
the equation 

dm’l = A{m^) X\dt), 

as an element o/Hq. Then its distribution converges locally weakly to the distribution of the 
solution path to the rough differential equation 

dmt = A{rnf) lK.{dt) 


started from xq . 


The interpolation result Theorem 2.1.2 follows directly from this fact, together with the 
construction of kinetic Brownian motion on AA as the image of its M'^-valued counterpart by 
Cartan’s development map, as described above, and the results of Section 2^ on the weak 

to the Stratonovich 


convergence of the rough path lift of kinetic Brownian motion on 
Brownian rough path. This closes the proof of the interpolation theorem 2.1.2 


The above result on transport of weak convergence only gives locally weak convergence of 
the image measures. The additional ingredient needed to turn that local weak convergence 
into genuine weak convergence, as stated in corollary 2.1.3, is a direct consequence of the 


following elementary lemma, in which we denote by Cr the set of continuous R/j-valued paths 
on [0,1] - recall Br stands for a geodesic ball of radius R. Denote by dA the boundary of a 
subset H of a topological space. 


Lemma 2.4.4. Let (Pn)n>o be a sequence of probability measures on Dq such that 

• the measures are all supported by C ([0, 1],AA), 

• they converge locally weakly to a probability measure P on Dq? which is also supported 
on C ([0, 1],A4.), and such that P(9Cij) = 0, for all R> 0. 

Then the sequence P„ converge weakly to F in C ([0, 1],AA). 
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Proof. First note that given any Borel set A of (^([0,1], A^) such that P(c}A) = 0, we have 

|P,(A) - P(A) I < |P,(A n Cr) - P(A n Cij) I + |P„(A n I + |p(a n c^) |. 

Since C([0,1], Al) = U_r>o^K) and since P is supported on (^([0,1], Al), given any £ > 0, we 
can hnd R large enough such that |l — P(Cij)| < e; in particular |P(y4 fl C^)| < e. We have 
also 

|Pn(AnC^)| < |i-p„(Cr)|. 

Moreover since P(5 Ch) = 0 we obtain that P(c}(A fl Cr)) = 0 and since P^ converges locally 
weakly to P we have 

p„(rp(>inc„)) ^ p(rs‘(-4nc„)). 

n —>+00 

Since T^^{A n Cr) = A D Cr it conies 

\Fn{AnCR)-FiAnCR)\ 0 

n ^+00 

and for the same reason P„(Cj?) converges to P(C'i?). Thus limsup^ |l —P„(C^) | < £ and hnally 

limsup |Pn(A) — P(v4)| < 2e 

n^oo 

which ends the proof of the lemma. □ 

3. Asymptotics on rotationally invariant manifolds 

We aim, in this last section, to try and understand the long time asymptotic behaviour of 
kinetic Brownian motion when t goes to inhnity and a is hxed. As for classical Brownian 
motion on a general Riemannian manifold, there is no hope to fully determine the asymptotic 
behaviour of kinetic Brownian motion on arbitrary Riemannian manifolds, as even for Brow¬ 
nian motion it is likely to depend on the base space in a very sensitive way; see for instance 
the work |ATU09j of Arnaudon, Thalmaier and Ulsamer for the study of the asymptotic be¬ 
haviour of Brownian motion on Cartan-Hadamard manifolds. Of crucial importance in the 
latter study is the fact that the distance to a hxed point dehnes a one-dimensional subdiffusion. 
The difficulties are actually a priori greater here as kinetic Brownian motion does not live on 
the base manifold, but on its unit tangent bundle, so that it is basically (2d — 1)—dimensional 
when the base manifold have dimension d, and there is no general reason why it should have 
lower-dimensional sub diffusions. 

The question of the long time asymptotic behaviour of a manifold-valued diffusion process is 
of course of different nature depending on whether or not the underlying manifold is compact. 
If it is compact, the question consists mainly in studying the trend to equilibrium of the 
process, and in relating the rate of convergence to the geometry of the manifold or to the 
parameters of the process. In our case, the inhnitesimal generator of kinetic Brownian motion 
is hypoelliptic and T^Ai is compact, if A1 is compact, and it is elementary to see that for all 
times t strictly greater than the diameter of Al, the density Pt{-,-) of the process is uniformly 
bounded below by a positive time-dependent constant. This ensures that the law of the process 
converges exponentially fast to the equilibrium measure, which is the Liouville measure on the 
unitary tangent bundle here. The question of determining the exact rate of convergence to 
equilibrium, in terms of the geometry of Al and the parameter a is difficult and will be the 
object of another work by the authors. 

If the base manifold is non-compact, the question of the long time asymptotics of a diffu¬ 
sion process amounts to hnd whether it is recurrent or transient, to exhibit some geometric 
asymptotic random variables associated to the sample paths, and, in the nicest situations, to 
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determine the Poisson boundary of the process. In order to get some tractable and signihcant 
information on the asymptotic behaviour of kinetic Brownian motion, we will restrict ourselves 
here on the case where the underlying manifold M is rotationally invariant. As in the case 
of classical Brownian motion, symmetries simplihes greatly the study as they allow to exhibit 
some lower dimensional subdiffusions of the initial process. The class of rotationally invariant 
manifolds is nevertheless rich enough to have some good idea of the interplay between geome¬ 
try, through curvature of the manifold, and the asymptotic behaviour of the process. Our goal 
here is to obtain some natural conditions on the curvature to ensure transience or recurrence 
of kinetic Brownian motion; if transient to determine its rate of escape, the existence or not 
of an asymptotic escape angle, and if possible to characterize the Poisson boundary of the 
process. We compare these conditions and results to those obtained for Brownian motion on 
rotationally invariant manifolds, in order to highlight the similarities and differences between 
the two processes. 


We recall in Section 3T the setting of rotationally invariant manifolds and present our results 
on the asymptotic behaviour of kinetic Brownian motion in Section 3^, comparing theorem 
to their Brownian analogue. As said above, we take advantage of this special symmetric 
setting to exhibit some subdiffusions that are easier to analyse than the full process. The 
radial component, together with its derivative, provide a 2-dimensional subdiffusion whose 
asymptotic behaviour is thouroughly studied in Section ^ As a result, we are able to decide 
on the recurrent/transient character of kinetic Brownian motion, and in the latter case, to 
exhibit a whole range of behaviours for its escape rate, depending on the expansion factor in 
the warped product defining the metric. The angular component is studied in Section |3.5 


while we show that the asymptotic angle generates the invariant sigma field of the process 
under rather general geometric assumptions. Working with a hypoelliptic diffusion makes that 
task highly non-trivial. 


3.1. Rotationally invariant manifolds. Let us recall that a rotationally invariant Rie- 
mannian manifold (A4,g) is a Riemannian manifold such that M. oi M. minus a point, admits 
a global polar coordinates system in which the metric g has a warped product structure. 
Namely, as illustrated in Figure below, (A4,g) is a rotationnaly invariant if M or Al\{o}, 
where o is a point in M, is diffeomorphic to the product (0, -|-cxd) x endowed with its 
standard polar coordinates (r, 0) in which the Riemannian metric g takes the form 

g = dr^ + f{rYd6'^, (r, 6) G (0, -|-cxd) x 


where / is a positive smooth function on (0,-1-cxd), say at least satisfying the classical 
assumptions /(O) = 0 and /'(O) = 1, see e.g. |Sto69j p. 179-183, and dO stands for the metric 
on inherited from the Euclidean metric of M'^. We refer to |GW79j for a detailed study 
of such manifolds. Since we are interested in non compact manifolds here, we shall suppose 
in the sequel that /(r) goes to infinity as r increases indefinitely. We shall assume in the 
sequel that d > 3; the case d = 2 presents no difficulty but requires a separate treatment, 
see Remark 3.3.3 in Section 3^ below, and Appendix A. 

The radial curvature K on a rotationally invariant manifold A4 is given by 


K{r) 


fjr) 

fir) ’ 


so that the curvature properties of the manifold are intimately related to the convexity of the 
warping function /. For example, global convexity of / is equivalent to global non positive 
curvature, and in that case we have f'ir)/f{r) > 1/r for all r > 0, in particular / is increasing. 
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Figure 4. Global polar coordinates system on a rotationally invariant manifold. 


Indeed, we have 


t 

f 





so that f”{r) > 0 



< 



2 


and integrating the ratio, since /(O) = 0 and /'(O) = 1, we get f'{r)/f{r) > 1/r. Moreover, 
if the warping function / is log—concave, the logarithmic derivative f /f is decreasing and 
converge to £ E [0, +cxd) when r goes to inhnity. In that case, for e > 0 and for r sufficiantly 
large, we have thus —{£ + < K{r) < 0. In particular, if / is log —concave and £ = 0, the 

rotationally invariant manifold is asymptotically flat. 


Example 3.1.1. Here are some classical examples of rotationnaly invariant manifolds and 
their associated radial curvature. 

(1) The d-dimensional Euclidean space minus a point is of course a rotationally invariant 
manifold. It corresponds to the choice f{r) = r for the warping function and has 
constant zero curvature. 

(2) The d-dimensional hyperbolic space EI'^ minus a point, seen as the pseudo-sphere in 
Minkowski space corresponds to the choice f{r) = sinh(r). The radial curvature 
is constant equal to —1. 

(3) If f is of polynomial growth, say f{r) = for r large enough, then the curvature 

K{r) = — has the sign of 1 — /3 et goes to zero as r goes to infinity, i.e. the 

manifold is asymptotically flat. 

(4) If f is of subexponential growth, say /(r) = exp (r^) with 0 < fl < 1 forr large enough, 
then the curvature 

flifl - 1 ) 

j.2-13 ^2(l-/3) 

is asymptotically negative and goes to zero as r goes to infinity. 

(5) If the growth of f is more than exponential, e.g. f{r) = exp (r^) with fl > 1 for r 
large enough, then radial curvature is negative and goes to minus infinity as r goes to 
infinity. 
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3.2. Statement of the results. Before expliciting the long time asymptotic behaviour of 
kinetic Brownian motion on M in terms of the geometry of M, let us recall the corresponding 
results for the classical Brownian motion (Xt)o<t<T on M. 


3.2.1. Asymptotics of classical Brownian motion. On a rotationally invariant manifold, the 
classical Brownian motion fortunately admits a one-dimensional diffusion, namely the radial 
subdiffusion (rt)o<t<r when written in polar coordinates (r^, Ot)Q<t<T- Moreover, the angular 
component {Ot)o<t<T can be shown to be a time-changed spherical Brownian motion on 
the clock being a simple additive functional of the radial process, see Example 3.3.3 of |Hsn 00]. 
The following classical results are thus simple consequences of the study of this radial one¬ 
dimensional diffusion, see for instance |attt] and the references therein. 

Theorem 3.2.1. Let {rt,9t)o<t<T stand for the polar decomposition of the classical Brown¬ 
ian motion with generator on a rotationally invariant manifold (M,g) parametrized by 

((0, -|-cxd) X dr^ p{r)d9^). 

(1) The lifetime r of the process is almost surely finite or infinite; it is finite almost surely 
if and only if 

/ + 00 / ^ + 00 \ 

f~\r) I^J f~^{p)dpj dr < -kcx). 

(2) The radial process rt is transient a.s. if and only if 

+ O0 

f^~‘^{r)dr < - 1 - 00 . 

(3) The angular process 9t converges a.s. to a random point 9^0 in if and only if 

/ +00 / I'+ca \ 

f-\p)dp\ dr <+ 00 . 

When the process is transient and admits a random escape direction 9^0 G it can 

be shown that this direction is the only non-trivial asymptotic variable, namely the Poisson 
boundary of the process coincides almost surely with ^(^oo), see e.g. Theorem 4 of jAT14] 
for a simple proof. Recall that the Poisson boundary of a Markov process can be defined 
either as its invariant sigma field or equivalently, via the classical bijection, as the set bounded 
harmonic functions of its infinitesimal generator. 

Theorem 3.2.2. Let {rt,9t)t>o stand for the polar decomposition of the classical Brown¬ 
ian motion with generator ^Aj^ on a rotationally invariant manifold (M,g) parametrized 
by ((0, -l-oo) X dr"^ -\- p{r)d9‘^) such that 

/ + 00 / ^+00 \ 

U f^-\p)dp\ dr < + 00 . 



Then, the Poisson boundary of the process {rt,9t)t>o coincides almost surely with cr(6'oo)- 

3.2.2. Asymptotics of kinetic Brownian motion. The following results sum up the long time 
asymtotics of kinetic Brownian motion on rotationally invariant manifolds. For simplicity 
and clarity of the statements, we give here the results under non-optimal assumptions. More 
precise and stronger results are given in Sections 3T, 3A and 3T, along with the proofs of the 
corresponding statements. 


As shown in Proposition 1.2.2 the lifetime of kinetic Brownian motion on any complete 
Riemannian manifold is almost surely infinite. The recurrent and transient character of kinetic 
Brownian motion are clarified by the following statement. 
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Theorem 3.2.3. Let {zt)t>o = {xt,Xt)t>o be the kinetic Brownian motion with values in 
the unitary tangent bundle T^Ai of a rotationally invariant manifold (M,g) parametrized by 
((0, +oo) X dr'^ + p{r)d9^) and let {xt)t>o = 9t)t>o be the polar decomposition of the 
first projection. 

(1) If K <0, then the process rt is transient almost surely. 

(2) If K >0 the process rt is transient almost surely if and only if 

+ 0O 

f^~'^{r)dr < +CX). 



Remark 3.2.4. Note that the condition ensuring the transience of the radial component of the 
kinetic Brownian motion is the same as the one for the classical Brownian motion. Indeed, the 
strategy of the proof of Theorem 3.2.3 , which is given in Section 3^ below, is to use a change 
of variables and a time change for the kinetic Brownian motion to make it look “similar” to 
its classical analogue and then use comparison results. 


The next proposition gives a sufficient condition for the convergence of the angular compo¬ 
nent of kinetic Brownian motion in terms of the warping function: for simplicity, we assume 
here that it is log —concave. We strongly believe that this condition is necessary but despite 
our efforts, we did not manage to give a proof of the corresponding statement, see Remark 
[3XT] below. 


Proposition 3.2.5. Let {zt)t>o = {xt,Xt)t>o stand for kinetic Brownian motion with values 
in the unitary tangent bundle T^M. of a rotationally invariant manifold {Ai,g), parametrized 
by ((0, -f-oo) X p{r)de‘^). Let {rt, 6t)t>o stand for the polar coordinates of {xt)t>o- 

Suppose that the warping function f is log—concave and satisfies 

/ + 00 / r+OO \ 

f-^{p)dpj dr <+^. 

Then, the angular process {9t)t>o converges almost surely to a random point 9oo on 


Remark 3.2.6. Note that this time, we obtain a convergence criterion which is different 
from the one concerning the classical Brownian motion. Namely, the exponent d — 3 of the 
classical case is replaced by d — 2 here. Therefore, if the integrability condition of Proposition 
3.2.5 is indeed necessary, there are situations where the angle 9t is convergent for the classical 
Brownian motion and does not converge for the kinetic Brownian motion. For example, if the 
warping function is of the form f{r) = r^ with (I > 0, then the angular component converges 
in the classical case if and only if (3 > 1, whereas it would converge in the kinetic case if and 
only if (3 > 2. 


Let us conclude the synthetic description of the long time asymptotic behaviour of kinetic 
Brownian motion by expliciting its Poisson boundary, as we did in the Brownian case, in 
Theorem 13.2.21 above. 


Theorem 3.2.7. Let {zt)t>o = {xt,Xt)t>o stand for kinetic Brownian motion with values in 
the unitary tangent bundle T^M. of a rotationally invariant manifold (N4,g), parametrized 
by ((0, -f-oo) X dr"^ + p{r)d9‘^). Let {rt, 9t)t>o stand for the polar coordinates of {xt)t>o- 
Suppose that the warping function f is log—concave and satisfies 

/ -l-oo / p+oo \ 

f-\r) p-^p)dpj dr < -Foo. 

Then the Poisson boundary of the whole process {zt)t>o coincides almost surely with a{9oo). 
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The proof of Theorem 3.2.7 which is given in Section 3T below, is based on the devissage 
method introdnced by two of the anthors in |AT14] . and which allows to identify Poisson 
bonndaries of diffnsion processes with valnes in manifolds provided the latter admit enough 
symmetries so that the original diffusion has some natural subdiffusions. See Section 4 of 
|AT14j for various examples of application of the method. This method is particularly well 
suited here since the warped product structure of the manifold gives automatically the exis¬ 
tence of lower dimensional subdiffusions; see Section 3.3 below. To help the reader by keeping 
the present paper as self-contained as possible, let us recall the main result of jAT14j . 


Theorem 3.2.8 (Devissage method - Theorem A of |AT14] b Let M he a differentiable 
manifold and G be a finite dimensional connected Lie group. Let {xt,gt)t>o be a diffusion 
process with values in J\f x G, starting from {x,g) G A/" x G. Denote by P(x,g) its law. Let us 
suppose that following conditions are satisfied: 

(1) Devissage condition: the process {xt)t>o is itself a diffusion process with values in Af. 
Its own invariant sigma field Inv{{xt)t>o) is either trivial or generated by a random 
variable ioo with values in a separable measure space {S,Q,X) and the law of ioo is 
absolutely continuous with respect to X. 

(2) Convergence condition: the second projection {gt)t>o converges almost surely to a ran¬ 
dom element goo ofG. 

(3) Eguivariance condition: the infinitesimal generator C of the diffusion {xt,gt)t>o is 
eguivariant under the action of G on G°°{J\f x G, M), i.e. V/ G G°°{Af x G, M), we 
have for {x, g, h) & AT x G x G: 

^f{-, 9 -)){x,h) = {Cf){x,gh). 

(4) Regularity condition: all bounded C—harmonic functions are continuous on the state 
space AT X G. 

Then, Inv{{xt, gt)t>o) o^nd Inv{{xt)t>o) V (T{goo) coincide up to ¥(^,^^g) — negligeable sets. 


Note that this result can be generalized to the case where the second projection does not take 
values in a hnite dimensional Lie group G but in a hnite dimensional connected co-compact 
homogeneous space Y := G/K, see Theorem B of |AT14] . 


3.3. Kinetic Brownian motion in polar coodinates. The remainder of this article is 


devoted to proving (hner versions of) Theorem 3.2.3, Proposition 3.2.5 and Theorem 3.2.7 


We shall exclusively be working, from now on, on a rotationnaly invariant manifold and use 
polar coordinates. In order to make explicit the long time behaviour of kinetic Brownian 
motion, we hrst write down the system of stochastic differential equations satished by the 
polar coordinates of kinetic Brownian motion on a rotationally invariant manifold. We exhibit 
in particular a lower dimensional subdiffusion that will facilitate its study. Let us denote by 
{r,6,f,6) the polar coordinates in T^AA. Then, system (1.3) giving the dynamics of kinetic 
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Brownian motion in any local chart reads here 

drt = hdi, 


(3.1) 


dOt = 9tdt, 


drt = adM[ - ^{d - l)rtdt + y (rj (l - r^) dt, 
^ J 




a 




dOl = adM^ - {d- l)eidt - 2^{rt)rtdldt - 61 

2 / 


Pin) 


dt, i = 1,... ,d, 


where and Mf are local martingales whose covariance matrix is given by 


(3.2) 


' d{M\M^)t ={1-Pt)dt, 

d{M\M^^)t =-hm, J = l, 

's, 


,d, 


d{M\M^')t = 




oiei 


Pin) 


6161] dt, i,i = l,...,d. 


Since the sample paths are parametrized by the arc-length, we have also the following metric 
relation, relating the radial and angnlar components: 

( 3 - 3 ) p + f{rtf\6t\^ = 1. 

Dne to isotropy, kinetic Brownian motion on a rotationally invariant manifold fortnnately 
admits a radial snbdiffnsion. Contrary to the case of the classical Brownian motion where the 
radial snbdiffnsion is one dimensional, the radial snbdiffnsion of kinetic Brownian motion is 
2 —dimensional. 

Proposition 3.3.1. The process {rt,ft)t>o with values in (0,-1-cx)) x [—1,1] is a subdijfusion 
of the kinetic Brownian motion {rt,ft,6t,6t)t>o in PM.- 


Proof. From Equation (3.1), the 2—dimensional radial process is clearly a subdiffusion of the 


whole process. Namely, there exists a real standard Brownian motion B such that {rt,ft) is 
solution to the following stochastic differential equations system: 

drt = hdt. 


(3.4) 


drt = a pi- pdBt - y (d - phdt + y (r^) (l - p) dt. 


□ 


Let us describe the behaviour of this radial process {rt,ft) at the boundary points of its 
state space (0, -|-oo) x [—1,1]. 

Proposition 3.3.2. If d >3, for all starting point (ro,'fo) G (0, +oo) x [—1,1], the stochastic 
differential equation system (3. 4) admits a unique strong solution, which is well defined for all 
t > 0, and such that rt > 0 and —1 < rt < 1 for all t > 0. 


Proof. First remark that if Fq < 1 and tq > 0, all the coefficients in (3.4) being smooth for 
(r, r) G (0, -|-cxd) x (—1,1), it admits a unique stro ng so lution up to explosion. Moreover, if 
f-g = 1, since r i—)■ pi — is 1/2—Holder, Equation (3.4) also admits a unique strong solution 
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up to explosion and fj<l for arbitrary small t > 0 almost surely. From the metric relation 


(3.3), we have < 1 so that rt can not go to inhnity in hnite time almost surely. Otherwise, 


if ^5 < 1 and ro > 0, straightforward ltd calculus shows that 


- rf) = /^(ro)(l - ro^) exp -aH + (d - 3)a^ 


rds + 2a 


vA 


zdB,, 


From this, we deduce that, almost surely, the left hand side can not vanish in hnite time, in 
particular > 0 and — 1 < r* < 1 for alH > 0. □ 

Remark 3.3.3. In the case d = 2, standard comparison arguments show that the endpoints 
— 1 and +1 can he reached by rt in finite time almost surely, and that these points are instan¬ 


taneously reflecting. From relation (3.3), the norm l^l vanishes at such hitting times, so the 


renormalized process 0/\6\ can not he considered. This is the only reason why the case d = 2 
reguires a separate treatment. Nevertheless, the study can be done in the same line as in the 
case d>3, see e.g. Appendix A for a complete treatment of the case where the base manifold 
is the hyperbolic plane. 


From Equation (3.3), Proposition 


3.3.2 


ensures that \9t\ > 0 for all t > 0 almost surely. 
It will be usefull in the sequel to consider the normalized angular process {0t,0t/\9t\)t>Q with 
values in Starting from Equation (3.1), a direct calculation shows that this process 


satishes the following system of stochastic differential equations, for 1 < i < d: 


(3.5) 



rt \ .. (T" , , dt 


dt - {d — 2) 


\0t\^ 


7^ + cdNlji, 


where the local martingales are given by 

1 


dNl := ^ I dMf - 

\9t\ V l^tl "ir \0t\ 


so that their covariance matrix reads 
(3.6) 


d{Ni,N}) = \ 0 - OX - 


dt 


\9t\ \9t\ J ^ - rT 

Moreover, for 1 < i < d, we have 

ll/dMtdMF-Af-S- 
\e,\ F / \e.\ \o,\ 

Equivalently, the inhnitesimal generator of the full process {rt,ft,9t,9t/\9t\)t>Q is given by 


(3.7) (M/, iv;> = 2- [ {dMt dMf A duf ) ] = 0. 


(3.8) 


C - C(r,r) + ( 


Vi — 
fir) 


^2 
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where £(r,r) is the infinitesimal generator of the subdiffusion {rt,ft)t>o, namely 
C(^r,r) = rdr + j{r) (l - r^) + y (-(d - l)rdr + (l - r^) dj) . 

Note that if cr = 0, the above infinitesimal generator is nothing but the generator of the 
geodesic flow on T^M.. 

Remark 3.3.4. Note that local martingale = {N^,..., Nf) can be represented by a standard 
Euclidean Brownian motion Wt = (hh/,..., W^) in the following way: 

dNt = dWt - OtiOt, dWt) - 7^ ( 7^, dWt 

\0t\ \|6't| 

In particular, dNt is orthogonal to both 9t and y. 

3.4. Asymptotics of the radial components. We study in this section the recurrence and 
transience properties of the radial subdiffusion. We give in particular the proof of Theorem 
13.2.31 stated in Section 13.2.11 



3.4.1. Transience/recurrence of the radial process rt- We first show that the under the single 
hypothesis that the radial curvature K{r) is globally non-positive, the first projection is 


almost surely transient, which is the first statement of Theorem 3.2.3 


Proposition 3.4.1. Let {rt,rt)t>o be the unigue strong solution of Eguation (S.f) starting 
from (ro,ro) € (0,-|-oo) x [—1,1]. If the warping function f is convex on (0,-|-oo) i.e. if the 
radial curvature K{r) is globally non positive, then r* goes almost surely to infinity with t. 

Remark 3.4.2. Let us stress that we are dealing here with global convexity and not convexity 
outside a compact set. Even in the case of the standard Brownian in a rotationally invariant 
manifold, negative curvature outside a compact set is not sufficient to ensure transience of 
process, see Remark 2.4 of 


Proof. Let {rt,rt) be the solution of Equation ( |3.4[ ) starting from (ro,ro). We have then 
2(ro - ft) ^ _ 


(3.9) 


n - ro = 



(d — l)cr2 


4(r,) (1 -r^) ds- 


{d — l)a 


^/l^sdBs 


:=Jt 


-V- 

:=Mt 


Note that the process R is bounded since ft is. Otherwise, by Itb’s formula, we have 

(3.10) 

rt rt j!i rt 


O = ^0 


{d - 1 )(tH + da"^ / {l-fl)ds+2 / —{rs)fs{l-f‘^^ds+2(T / r^-^l - fldBg 
Jo Jo / Jo 


= {M)t 


—V— 

■■=Kt 


—V— 

:=Nt 


Let us first remark that the martingale Mt and the process Jt can not converge simultaneously 
when t goes to infinity. Indeed, if it was th e cas e, since f^ < 1, the processes {M)t, Kt and 
Nt would also converge and from Equation (3.10), f^ would tend to minus infinity, hence the 


contradiction. Thus, if we suppose that Mt is convergent, necessarily Jt is divergent and by 

ipose now that Mt is not convergent so that the following 
^t are well defined for all f > 0 and both go to infinity when 


Equation (3.9), so is rt. Let us suppose now that Mt is not convergent so that the following 
time change Dt and its inverse 
t goes to infinity: 

Dt := a‘^{M)t = [ (1 - fl)ds, 


^ := inf{s > 0, > t}. 
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To simplify the expressions, let us define pt := Pt introduce the process 

{d — l)a^ 


Ut (jpt + Pt, 

Then, Ut satisfies the equation 

(3.11) 


where a := 


1 f 

+ dBt, 

f 


i.e. 


1 f' _ 

dut= — - Pt)) dt + dBt 


Now consider the process Vt starting from vq = mq and solution of the following stochastic 
differential equation 

, d-I dt , ^ 
dVs = -h dBt- 

2 1 + Uf 

Note that Ut + 1 is then a standard Bessel process of dimension d > 3, in particular it is 
almost surely transient. Otherwise, as already noticed at the end of Section 3T, thanks to the 
curvature hypothesis, we have f /f{r) > 1/r. Recalling that pt G [—1,1], we deduce that 


d{vt - ut)^ = 


d-l 


< 


2 

d-l 

2 

d-l 


(?-■(«,-A)) 


^Vt>Ut 


dt 


1 


l+Vt Ut-pt 

{ut - Vt) - (1 + pt) 
{I + Vt){ut - pt) 


Vt>Ut 


dt 


1 


dt < 0. 




In other words, we have Ut > Vt almost surely. We can then conclude that Ut is transient, and 
since pt is bounded, pt is transient and finally the time changed process = p^^ is transient. 

□ 

Remark 3.4.3. Ifrt is transient almost surely and if moreover the warping function f is such 
that f'{r)/f{r) goes to zero when r goes to infinity, then Eguation (3.10) shows that almost 
surely 


1 r 

hm T / (1 
t —>-+00 i Jq 


fi)d, = 1 - 


i.e. 


hm - 

t —^“1“00 ^ 




'0 


As in the case of Brownian motion on rotationally symmetric manifolds, the transience 
of rt can also be caracterized by the integrability of the inverse of the function /. Since 
the radial subdiffusion is 2—dimensional here, things are less simple. Nevertheless, we can 
charaterize transience of kinetic Brownian motion under a simple monotonicity condition on 
the logarithmic derivative of /. 


Proposition 3.4.4. Let {rt,rt)t>o be the solution of Eguation (S.f) associated with the func¬ 
tion f and starting from (ro,ro) G (0,+oo) x [—1,1]. 

(1) Suppose that there exists a smooth function g on (0,+cx)), such that g is positive and 
log—concave, and such that f'/f > g'/g. Then, if g^~^ < +cxd, the radial process 
rt is transient almost surely. 

(2) On the contrary, suppose that there exists a smooth function h on (0,+cx)), such that 
h is positive and log—concave, and such that h'/h > f'/f. Then, if h^~^ = +cx), 
almost surely the radial process rt is not transient. 

(3) In particular, if the warping function f is log—concave, then the radial process rt is 
almost surely transient if and only if f^~'^ < +cx). 
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Proof. Let {rt,ft) be as in the above statement. With the same notations as in the proof of 
Proposition 3.4.1, the process Ut is then solution of Equation (3.11), namely 

^ f ~ 

dut = —— - pt)) dt + dBf 




f 


Let us hrst suppose that there exists a smooth function g on (0, +oo), such that g is positive, 
log—concave, and such that f'/f > g'/g. Consider the process Vt starting from vq = Uq and 
solution of the stochastic differential equation 

1 0^ 

dvt = ^— Cd~^{vt + 1 )) dt + dBt. 

<jg 

Classical comparison results then show that Ut > Vt almost surely for all f > 0. To conclude, 
note that the process Vt := + 1) is solution of the equation 


dvt = {vt) dt + dBt, 

2 g 

and Vt is transient if (and only if) Ji°° g^~^ < +C) 0 , see e.g. Theorem 1.1 p. 208 of |Pin95j . 
hence the first point. Suppose now that there exists a smooth function h on (0,+cx)), such 
that h is positive, log —concave, and h'/h > f'/f. As above, classical comparison results then 
show that Ut < Vt almost surely for all f > 0, where Vt is now the solution starting from uq = '^o 
of the equation 

1 h' _ 

dvt = {o-~^{vt - 1)) dt + dBt. 

Again, to conclude, remark that new rescaled process Vt := — 1) is bow solution of 

the equation 

7~ d—lh' 

dvt = ^ — [Vt) dt + dBt, 

and Vt is recurrent if (and only if) = +cxd, hence the result. Now if the warping 

function is log —concave, combining the first two points, we deduce that r+ is transient if and 
only if = +CX). □ 


The last result allows to caracterize the transience of the radial process in non negatively 
curved rotationally symmetric manifolds. Indeed, recall that iP = —f"/f = —[f'/fY—if'/f)'^, 
thus if iP > 0, we have necessarily [f'/f)' < 0 i.e. / is log—concave. From Proposition 3.4.4 


we thus deduce the following results which is the second point in Theorem 3.2.3 


Corollary 3.4.5. If the radial curvature K is globally non negative, then the radial process rt 
is transient if and only if < +cxd. 

3.4.2. Transience/recurrence of rt. We now describe the asymptotic behaviour of the radial 
derivative rt in the case where rt is transient. We first consider the simplest case where the 
ration f'/f is constant and thus rt is a one dimensional diffusion in [—1, +1]. We then have 
the following lemma: 


Lemma 3.4.6. Let [rt,rt)t>o be the solution of Equation (d.f) associated with the warping 
function f such that f / f[r) = £ G M and starting from (ro, ro) G (0, +oo) x [—1,1]. Then the 
process ft is ergodic in the interval (—1,+1) with invariant probability measure 


Pi[dx) := 


[1-x^) 


d-3 
2 


2^ . 


-dx. 


(1 


X 


d-3 

I 2 eo- 


21 . 


dx 


'-1 
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In particular, if £ > 0 we have almost surely 

n 


lim — = / Xiae{dx) G (0, +C)o). 

t —i j ^ 


Proof. If f'/f is constant, then Equation (3.4) shows that is a one-dimensional diffusion 
and the probability measure whose support is the interval (—1, +1), is invariant, hence the 
result. The ergodic theorem then ensures that r^ is ballistic almost surely. □ 

Now if f'/f is not constant but bounded outside a compact set, we can then deduce that 
Vt is Harris recurrent as soon as r* is transient. 


Proposition 3.4.7. Let {rt,rt)t>o the unique strong solution of Equation (S.f) starting 
from (ro,ro) G (0, -1-cx)) x [—1,1]. Suppose that there exists R > 0 and £ > 0 such that for all 
r > R, we have \f'{r)/f{r)\ < £. Then, ifr^ is transient almost surely, the process ft is Harris 
recurrent in (—1,-|-1). 

Proof. If rt is transient almost surely, using classical one dimensional comparison results, we 
have ef < ft < ef for t sufficiently large where the two processes ef are solutions of the 
equations 


def = a 


fdBt 


a 


— {d— l)efdt ± £ (l 


.±|2 


) dt. 


From Lemma 3.4.6, both processes are ergodic in (—1, -fl), hence the result. 


□ 


Corollary 3.4.8. Let {rt,ft)t>o be the unique strong solution of Equation (3. 4) starting 
from (ro,ro) G (0,-|-oo) x [—1,1]. Suppose that the warping function f is log—concave with 
hm,._,.+oo f /fiff") = £ > 0. Then the process ft is Harris recurrent in (—1, -f-l). 

Proof. Using the monotonicity of f'/f, we have f'/f{rt) > £ almost surely. By classical 
comparison theorems, we deduce that ft is bounded below by its analogue with f'/f = £. In 


particular, by Lemma 3.4.6 , we obtain that rt is transient almost surely. Then using again 
comparison results and Lemma 3.4.6 we deduce that ft is Harris recurrent. □ 


The almost sure asymptotic behaviour of ft is not clear in the case where f / f is not bounded 
above but does not go to inhnity. Nevertheless, if the logarithmic derivative is non decreasing 
outside a compact set and goes to inhnity with r, then one can show that ft converges to one in 
probability. A typical example of this last case is a smooth function / such that /(r) = exp(r^) 
for r > 1 with /3 > 1 so that f /f{r) = grows to inhnity with r. 

Proposition 3.4.9. If f is convex and if f /f is non decreasing outside a compact set and 
goes to infinity with r, then ft converges in probability to one as t goes to infinity. 


Proof. If / is convex, we know that rt is transient almost surely by Proposition 3.4.1 In 


particular, since f /f{r) goes to inhnity with r, by classical comparison results, we get that 
for all £ > 0 and for all t sufficiently large, ft > where is the solution of Equation 


(3.4) associated with a warping function whose logarithmic derivative is constant equal to £. 


Therefore, if £ > 0, for t sufficiently large, we have 
F{ft > I — e) > F[uI > 1 — e) 


t —^ ~1“CXD 




where the last convergence is a consequence of Lemma 3.4.6| and of the ergodicity of the process 
(mj). a direct calculation shows that — e, 1]) goes to one as £ goes to inhnity, hence the 
result. □ 
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3.4.3. Rate of escape. We conclude this section dedicated to the study of the radial process 
by expliciting the rate of escape of r* when transient. Under some simple extra assumptions 
on the behaviour of /'// at infinity, it is indeed possible to caracterize the speed of divergence 


of Tt. With the same notations as in the proof of Proposition 3.4.1, let {rt,rt) be the solution 
of Equation (3.4) associated with the function / and starting from (ro,fo) G (0, +oo) x [—1,1] 
and Ut the associated solution of Equation (3.11). If we assume that / is log—concave, then 


the same arguments than the ones used in the proof of Proposition 3.4.4 show that the speed 
of divergence of a~^u^ 2 f- and thus the speed of divergence of is the same as the one of the 
process Vt solution of 


(3.12) 


cRji = 


d-l f 
2 / 


(vt) dt + dBt 


Recall that pt is obtained from r* by the simple time change 

Dt = a^ [ (1 - rl)ds, 

Jo 


and from Remark 3.4.3 we know that almost surely, the clock Dt is asymptotically linear in 


t as soon as f /f{r) goes to zero at inhnity. In hne, if f /f{r) decrease to zero as r goes to 
inhnity, up to a scalar factor, the rate of escape of rt is the same as the one of Vt solution 
of Equation (3.12) given above. So kinetic Brownian motion happens to escape to oo at the 


same speed as Brownian motion, up to a multiplicative constant. 


Proposition 3.4.10. Let {rt,ft)t>o be the unique strong solution of Equation (3.f} starting 


from (ro,ro) G (0,+oo) x [—1,1]. In the following cases, up to a multiplicative constant, the 
rate of escape of the process rt is the same as the one of the solution Vt of Equation (3.12) .- 

(1) If f is of polynomial growth at infinity, namely if f{r) = r^ for r large enough with 
{d—l)f > 1, then the process Vt behaves asymptotically as a Bessel process of dimension 
d' := 1 + (3{d — 1). 

(2) If f is of subexponential growth at infinity, e.g. if f{r) = exp(r^) for r large enough 

^ 1 

with 0 < /3 < 1, then Vt goes to infinity at speed t'^-f^, in the sense that the ratio of the 
two terms converges almost surely to a positive deterministic explicit constant. 

(3) If f is of exponential growth at infinity, e.g. /(r) = for c> 0, then rt/t converges 
almost surely to a positive deterministic explicit constant. 


Proof. In the three cases, we know from Theorem 3.2.1 that the process Vt is almost surely 
transient. In the hrst case where /(r) = r^. Equation (3.12) reads 


{d-l)fidt 

dwt — ---— + dBt, 

2 Vt 

i.e. Vt is a Bessel process of dimension d' = {d — 1)(3 + 1, hence the result. In the second case 
where /(r) = exp(r^) for r large enough, we have for t large enough 

(3.13) dvt — - dBt, 

2 Vt ^ 

or equivalently if H{x) := : 

(3.14) dHifvt) = (^{d - 1)13 + dt + (2 - (3)vl~^dBt 
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Since Vt is transient almost snrely, from Eqnation (|3.14|) we have 
(3.15) 


H{v,) = I + o(l) I t + M„ 


where Mt is a martingale with bracket 

{M)t = (2 




r2- 


^^ds. 


Otherwise, from Eqnation (3.13), we have that Vt = o(t) almost snrely as t goes to infinity. 
Injecting this first estimate in the expression for the bracket of Mt, we get that Mt/t converges 
almost snrely to zero as soon asl/2</3<l, in other words, as t goes to infinity, we have 
almost snrely 


(3.16) 


, 2-/3 


lim 

t^ + OO 


(3{2-f3)id-l) 


Now if 0 < /S' < 1/2 < /3, classical comparison resnlts ensnre that the solntions Vt and v\ of 
(3.13) associated with fd and fd' respectively, and starting from a same point, satisfy v't < Vt 
almost snrely for all f > 0. As before, we have 


(3.17) 


where 


H{v't) = 


/3'(2-/S')(d-l) 


+ o(l) ) t + 




/\2 


t2- 


^^'ds. 


Since (d > 1/2, we know that almost surely Vt = and injecting this new estimate 

in the expression of {M')t, we get that M^/t goes almost surely to zero as t goes to infinity 
as soon as 1/4 < fd' < 1/2. Iterating the above argument, we eventually obtain that the 
almost sure convergence (3.16) holds true for all 0 < ,5 < 1. Finally, in the third case where 


the warping function has exponential growth, the logarithmic derivative in Equation (3.12) is 
constant, and Vt is then trivially ballistic, hence the result. □ 


3.5. Asymptotics of the angular components. Let us now describe the asymptotic be¬ 
havior of the angular components in function of the geometry of the underlying manifold. 
From Equation (3.3), we have 


0t = 00 + 


Thus, if we introduce the new clock Ct : = 


Os 


j,2 




-ds. 


\/l 

f{rs) 


-ds with inverse C* we have 


0t = 00 + / • “ ds. 

Since the last integrand is bounded by one, we deduce that 0t converges to a random point 0^0 
on the sphere as soon as Ct converges when t goes to inhnity. Therefore, a sufficient condition 
for the almost sure convergence of the angle is the almost sure hniteness of Ct as t goes to 
infinity. 
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Remark 3.5.1. It is tempting to say that the finiteness of Ct as t goes to infinity is also a 
necessary condition for the almost sure convergence of the angle. For example, if Ct goes to 
infinity and 9t converges to 6 qo, in Eguation ( |3.5 ) describing the evolution of 6 t/\dt\, the first 
term of the right hand side would go to infinity in the direction of O^o whereas the last terms 
would be asymptotically orthogonal to 6^0 (see Remark 3.3.4), which seems in contradiction 
with the boundedness of the left hand side. Unfortunately, we were not able to give a clear 
proof of this fact. Nevertheless, the next proposition, which completes the proof of Proposition 
3.2.5 , gives a necessary and sufficient condition, in terms of the integrability of f, for the 
almost sure convergence of the clock Ct. 

Proposition 3.5.2. Suppose that the warping function f is log—concave and satisfies the 
integrability criterion f^~‘^{p)dp < +cxd, then the clock Ct converges almost surely if and 
only if 


r‘+oo 


W) ■■= 




r‘+oo 


/ 


l-d 


{p)dp^ dr < +CX). 


Proof. Recall that, by Proposition 3.4.4, rt is transient almost surely if “ is integrable. 


By hypothesis, the warping function /(r) goes to inhnity with r. If it is log—concave, the 
logarithmic derivative f /f{r) then admits a limit £ > 0. Let us hrst suppose that £ > 0. In 


that case, / has exponential growth and classical comparison results and Lemma 3.4.6 show 
that rt is ballistic. Therefore, the process Ut converges almost surely as t goes to inhnity and 
the integral Id{f) is of course convergent, hence the result. Suppose now that i = 0. Applying 
Ito’s formula, we have 


d 


ViZ. 

fie) 


\/l - f 

fie) 


2 nf ( 
a - 2 j{rt)rt 


a‘^ f id — 2 )dt 

dt+^ — - 7= 

2 \/(rt)^r^ 


^.2 


— a 


e 


fie) 


dBt, 


and 


d 


V fie) J V fie) ) 


-da'^ - 3j{rt)rt 


, a 
dt + j 


[d — l)dt 

fie) 


-2a 


"t\/r 


fie) 


-dBt, 


Since rt is transient and f /fir) goes to zero as r goes to inhnity, from these two equations, 
we deduce the almost sure equivalences 


r*H-oo 


'*+00 


fie) 


■dt < +CX) 




'>+00 


fie) 


-dt < +CXD 


dt 


fie)V^ 


--dt < +CXD. 


In other words, with the notations of the proof of Proposition 3.4.1, the process Ct converges 
almost surely as t goes to inhnity if and only if 


'‘+00 


A 


'*+00 


fie) 


-dt = 


dt 

fipt) 


< +00 a.s., 


3.4.3 


the clock Dt 


where pt = and Dt = a^ f^(l — rg)ds. Otherwise, from Remark 
is almost surely asymptotically linear in t, and from the proofs of Propositions 3.4.1 and 
3.4.4 the process pt is asymptotically equivalent to the process a~^Ut, where Ut is solution of 


Equation (3.11). In hne, we have the almost sure equivalence 

/* + O0 


dt 


^+00 


fipt 


< +CX) 


dt 

fiVt) 


< +00, 
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where yt is the solution of the equation 

dyt = r, . at + dBt. 

2 f{yt) 

The almost sure finiteness of the last integral is then equivalent to the fact that the new time 
changed diffusion yt solution of 

dyt = ^^-^fiyt)dt + V f{yt)dBt, 

has an almost surely finite lifetime, which is finally equivalent to the finitness of Id{f), see e.g. 
Theorem 1.5 p. 212 of |Pin95j . □ 

3.6. Poisson boundary of the process. To conclude this section, we make explicit the 
Poisson boundary of kinetic Brownian motion in rotationally invariant manifolds whose warp¬ 
ing function / is log —concave and for which is integrable at infinity. We show that 
the invariant sigma field of the whole kinetic Brownian motion {zt)t>o = coincides 

almost surely with the sigma field generated by the escape angle when convergent. As an¬ 
nounced above, we take advantage here of the powerfull devissage method introduced in |AT14j 
that allows to compute the Poisson bourdary of a diffusion process starting from the one of 
a subdiffusion. Recall Theorem 3.2.8| giving its main range of application. This method is 
particularly well suited here because, thanks to the warped product structure of rotationally 
invariant manifolds, the radial process is always a subdiffusion of dimension two of the whole 
process. The plan of the proof of Theorem 3.2.7 is thus the following 


using coupling argu¬ 
ments, prove hrst that the Poisson boundary of the radial process (rt, ft)t>o is trivial, then use 
the devissage method to deduce that the Poisson boundary of the full process is generated by 
the escape angle. 


3.6.1. Poisson boundary of the radial subdiffusion. In this first section, we show that if 
is transient and is Harris recurrent, then the Poisson boundary of the radial subdiffusion 
is trivial. Let us first make a very simple but crucial remark: the first projection 
increases if and only if A is non negative so that the typical behaviour a of sample path of the 
radial subdiffusion looks like Fig. below. 



Figure 5. Typical behaviour of the sample paths of the radial subdiffusion. 

In this situation, we can exhibit a shift coupling between two independant copies of the 
radial process {rt,rt)t>o- Recall that the existence of shift couplings between sample paths 
is equivalent to the fact that the underlying process has trivial Poisson boundary, see e.g. 

|(iwnn] . 

Proposition 3.6.1. If rt is transient almost surely and if rt is Harris recurrent in (—1,-1-1), 
then the Poisson boundary of the radial subdiffusion {rt,rt) is trivial. 
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Proof. Let and be two independant solntions of Eqnation (3.4) starting from 

(rQ,ro) and (rQ,fg) respectively. Withont lost of generality, we can snppose that Tq < Tq. Let 
ns define S := inf{f > 0, rj = 1/2} and T := inf{f > S, rj = 0}. Since r/ is Harris recnrrent, 
S and T are finite almost snrely. Since rj increases if and only if rj is non negative, we have 
almost surely 

inf a{r) > 0, where a{r) := sup{r/, r/ = r} > 0. 


^'0 


r>r. 


t>0 


Geometrically, this means that the top of the grey zone associated with (r/,r/) on Fig. 
below is always non negative. Since r/ goes almost surely to infinity with t, the hitting time 
T' := inf{t > 0, r/ = rf} is finite almost surely. If > (y{rf), since r/ is recurrent, the 
curves (r/,r/) and (r/,r/) must intersect in finite time. Indeed, r/ must visit zero after T' and 
the grey zone is a barrier, this situation is illustrated in red on Fig. @ below. If rf < (y{r}p), 
then the two curves must have intersected before T, this situation is illustrated in green on 
Fig. § □ 



Figure 6. 


Typical behaviour of the sample paths of the radial subdiffusion. 


Corollary 3.6.2. Suppose that the warping function f is log—concave and satisfies the inte- 
grability criterion f^~^{p)dp < +cx), then the Poisson boundary of the radial subdiffusion 
{rtfft) is trivial. 


Proof. From Proposition 3.4.4 if is integrable at infinity, the process rt is transient almost 
surely. Otherwise, if / is log —concave, then /'// is bounded ouside a compact set and from 
Proposition 3.4.7, rt is then Harris recurrent in (—1, +1). Thus the assumptions of Proposition 
3.6.1|are fulhlled, hence the result. □ 


3.6.2. A diffusion lifting kinetic Brownian motion. Let us first observe that one cannot directly 
use the devissage method to make explicit the Poisson boundary of the kinetic Brownian 
motion {rt,rt,9t,9t/\0t\)t>o. Indeed, even if 9t converges almost surely to 6*oo when t goes 
infinity, the dynamics of 9t/\9t\ involves terms in 9t (to keep 9t on the sphere and the 

process (r*, ft, 9t/\9t\)t>o is not a subdiffusion although we do not expect invariant information 


from this term. In particular, the devissage condition of Theorem |3. 2. 8| is not fulhlled. In order 
to perform the setting of application of the devissage method we thus need to represent the 
angular process {9t/\9t\,9t)t>o € by a process with values in a bigger space where the 

devissage method can indeed be applied. Namely, in two steps, we will represent the angular 
process as a projection of some process {bt,gt)t>o living in SO{d — 1) x SO{d). Recall that 
(ei,... ,ed) denotes the canonical basis of and let us denote by SO{d — l) (resp. SO{d — 2)) 
the subgroup of SO{d) made of elements hxing ei (resp. hxing both ei and 62). 
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First remark that the angular process takes values in which can be 

seen as the homogeneous space SO{d)/SO{d — 2). Then consider the process {rt,rt,et)t>o, 
where e* = (e^,..., ef) G SO{d), whose inhnitesimal generator is given by 


(3.18) 

where 


Q = C^r,r) + a{r,r)Ho + y/?(r) ^ Id/, 

i=3 


a r,r : = 


\/l - 


f3{r) := 


j ej 62). In other words, the process {et)t>o is 


/(r) ’ 1 - r2 

and Ho^g) := gH^ (resp. Vj{g) = gVj) is the left invariant vector held in SO{d) generated by 
Hq := ei 0 62 - 62 ® Cl (resp. Vj := 62 0 e 
solution of the Stratonovich stochastic differential equation 

d 

(3.19) det = a{rt,ft)etHodt + f3{rt)'^etVj o dBl, 

i=3 


for a standard {d — 2)-Euclidean Brownian motion Doing so, the process 

(rt, rt, et(ei), et(e2))t>o has the same law as the target process (^rt,rt,dt,dt/\Ot\)^yQ i-e. in this 
hrst step, we lifted kinetic Brownian motion to a diffusion on M’*' x [—1,1] x SO{d). Let us 
now consider the time changed Brownian motion (&t)t>o on SO{d — 1) solution of 


d 

(3.20) dbt = /3{rt)Y,btViodBi. 

i=3 

Finally, viewing bt as the element diag(l,6i) of SO{d), dehne the process {gt)t>o on SO{d) by 
gt := etb^^ so that we have 

(3.21) dgt = a{rt, rt)gt h Hq b^^dt. 

Note that, by construction, we have gt{^i) = Ot and gtbt{e 2 ) = 0t/\9t\, i.e. after this second 
step, kinetic Brownian motion is now represented by a diffusion (r^, r*, bt, gt)t>o with values in 
M’*' X [—1,1] X SO{d — 1) X SO{d). But, now we are in good position to apply the devissage 
scheme. Indeed, the process {rt,ft,bt)t>o is a now subdiffusion of (j’t,ft,bt,gt)t>o- Moreover 
one shows that this diffusion has a trivial Poisson boundary. 


Proposition 3.6.3. Suppose that the warping function f is log —concave and satisfies the inte- 
grability criterion f^~‘^{p)dp < +00, the Poisson boundary of the subdiffusion {rt,rt,bt)t>o 
is trivial. 


Proof. Let us hrst remark that, following |CW00j . it is sufficient to construct a shift coupling 
between any two (r^, r^, fe^j-dihusions started from two diherent initial conditions to get the 
triviality of the Poisson boundary. We can construct such a shift coupling in the follow¬ 
ing way. Let (r^, and (r^, f^, be two independant versions of the [rt,rt,bt)- 

dihusion starting from {rl,rl,bl^ and respectively. We know from Proposition 


Py ,^2 ;,2)-almost surely, and such that one can modify the process [rj:, rj:, so as to have 

'f/i+t) = 'r/a+t)! 1 > 0. In particular, for all t > 0, we have 

ds _ ds 

Jt, Jt2 


3.6.1 


that there exists random times Ti,T2 that are hnite P(^i ,^1 (,i)-almost surely, respectively 
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Since Brownian motion in SO{d — 1) is ergodic, one can find two Brownian motions {c^{t))t>o 
and (c^(t))i>o on SO{d — 1), started from 6^^ and respectively, which are independant of 
the radial snbdiffusions and which conple almost surely in hnite time. Then, the processes yl, 
i = 1,2, dehned by the formulas 


yl 


rl,rl,bl), 


for 0 < f < Tj, 



for t > Ti, 


have the laws of (rt,-ft, -diffusions started from {rl,rl,bl) and (rQ,fQ, 
will couple in hnite time almost surely, hence the result. 


, respectively and 
□ 


3.6.3. Poisson boundary of the full diffusion. We can now give the proof of Theorem 3.2.7 
describing the Poisson boundary of the full kinetic Brownian motion. As mentioned before, 
the key tool used in this proof is the devissage method introduced in |AT14j . with the twist 
that we shall not apply it to kinetic Brownian motion, but rather to the lifted diffusion 
{rt, ft, gt, bt)t>o with values M+ x [—1,1] x SO{d) x SO{d — l) introduced in the last paragraph. 
In a hnal step one obtains the Poisson boundary of the kinetic Brownian motion from the 
Poisson boundary of the this lift. 


We have shown in Proposition |3.6.3] that if / is log—concave and satishes the integrability 
f^-‘^{p)dp < - 1 - 00 , then the Poisson boundary of the subdiffusion {rt,ft,bt)t>o 

is 


criterion 

is trivial. Moreover, if the second integrability condition Id(f) < oo of Proposition 3.5.2 
satished, the clock Ct is almost surely convergent and from Equation (3.21), the component gt 
then converges almost surely to poo in SO{d). In other words, the devissage and convergence 
conditions of Theorem 3.2.8 are satihed, with A/" = M+ x [—1,1] x SO{d — 1), G = SO{d), 
with (rt, ft, bt)^^Q in the role of the subdiffusion. Moreover, by equations (3.20) and (3.21), the 


inhnitesimal generator Qt, of the full diffusion {rt,ft,gt,bt)t>o is very similar to the generator 
Q dehned in (3.18), except the geodesic how is now conjugated, namely 


Gb ■= G(^r,r) + f{f)'^ yf + a(A f)Hb, 


i=3 


where Hb denotes the left invariant vector held on SO{d) dehned at point g by the formula 
gbHob~^ G TgSO{d). It follows that the equivariance condition (3) of theorem 3.2.8 is fulhlled, 
while it is easy to check that Qb satishes Hormander condition and, so, is hypoelliptic. As 


the regularity condition is also fulhlled, the devissage method applies. Namely, Theorem 3.2.8 
implies the following statement: 


Proposition 3.6.4. Suppose that the warping function f is \og-concave and satisfies 

/ +00 / ^- 1-00 \ 

f-^p)dpj dr < + 00 . 

Then, the Poisson boundary of the diffusion {rt, ft, h, gt)t>o coincides almost surely with O'{goo)- 

To hnally deduce the Poisson boundary of kinetic Brownian motion starting from the Poisson 
boundary of its lift {rt,ft,bt,gt)t>o with values in M+ x [—1,1] x SO{d) x SO{d — 1), it is 
convenient to introduce the map y from M’*' x [—1,1] x SO{d) x SO{d — 1) to the base space 
X [—1,1] X given by the formula 

x(r,f,g,b) = {r,f,g{ei),gb{e 2 )). 
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Let now / : x [—1,1] x —)■ M be a bounded harmonic function for the infinitesimal 

generator of kinetic Brownian motion {rt,rt,dt,Ot/\Gt\)t>o- Then, f o x is ^b-harmonic, and 
by Proposition 3.6.4 there exists a bounded measurable real-valued function F on SO{d) such 
that we have 

(/ ° X) {r, r, g, b) = E^r,r,g,b) [F(^oo)] , 
for all {r,f,g,b). Let now £x (^r,f,9,9/\9\) and g G SO{d), such that 

{r,r,9,9/\9\) = x{r,r,g,ld), 


and remark that we have = x(r, r, 5 ffc, fc for any k G SO{d — 1), so one can 

write 


(3.22) /(r, r, 9, 9/\9\) = E(^r,r, 9 k,k-^) [F(5foo)], 

for all k G SO{d — 1). The conclusion will then come from the following equivariance result. 

Lemma 3.6.5. The law of {gt)t>o under P(r,r,c,fc,fc-i) coincides with the law of {gtk\yQ under 
P(r,r,c/,id), for all k G SO{d — 1). 


Proof. Consider the processes {gt,bt)t>o starting at {gk,k~^), and set g[ := gtk and b'^ := kbt. 
Then, using equations (3.20) and (3.21), one easily checks that the process b[ starts from the 
identity and satishes the equation 


1=3 

while the process {g't)t>o starts from g and satisfies the equation 

dg't = a{rt,ft)g[b[HQ{b[)~^dt. 


So the law of the process (r^, rt, g[, giving the result. 


□ 


It follows from the previous lemma that the law of g^o under P{r,r,gfc,A:-i) coincides with the 
law of gook under P(r,r,g,id). However, as 


(3.23) 


/(r, r, 9, 9/\9\) = E(r,As,id) [F(^ooA;)] , 


for all k G SO{d — 1), from equation ( |3.22[ ), we have 
(3.24) 


f{r,r,9,9/\9\) = E(r,qs,id) 


P{gook) dk 


'SO{d-l) 


where dk stands for the normalized Haar measure on SO{d — 1). Now, given any measurable 
section S : ~ SO{d)/SO{d — 1) —>■ SO{d) to the natural projection, we dehne a bounded 

measurable real-valued function F setting 

F(0) := [ F{Si9)k) dk. 

JsO{d-l) 

In those terms, we hnally have 


f{r,f,9,9/\9\) = Ep,,.,g,id)[F(^oo(ei)) 


E 


(r,r,0/|e|,6)) 



which shows indeed that the Poisson boundary of kinetic Brownian motion is generated by 
the escape angle 9oo- 
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Appendix A. Kinetic Brownian motion in the hyperbolic plane 

In the preceeding section, we exhibited the almost sure long time asymptotic behavior of 
the kinetic Brow nian motion in rotationally invariant manifolds of dimension d > 3. As noted 
the case d = 2 requires a separate treatment since in that case the norm 


3.3.3 


in Remark 

vanishes almost surely in hnite time so that the process 9t/\9t \ is not well dehned. Nevertheless, 
this is just a formal difficulty due to the choice of polar coordinates, and an exhaustive study 
can be done in the same lines as above on a rotationally invariant manifold of dimension 2. 
To emphasize this fact, let us consider the case of the hyperbolic space of dimension 2. 

A.l. Kinetic Brownian motion in We consider the half-space model of the hyper¬ 

bolic plane, namely we view as the half-space {(x, y) G ?/ > 0} endowed with the metric 
ds^ = y~^{dx‘^ + dy'^). In this coordinates system, the kinetic Brownian motion {xt,yt,Xt,yt) 
started at {xo,yo,Xo,yo) G is the solution of the system of stochastic differential equa¬ 


tions (1.3) which simply writes here: 


(A.l) 


dxt 

= Xtdt, 



dyt 

= ytdt, 



dxt 

- 

cr^ . 

—Xtdt 

a dMf, 

yt 

2 


dyt 

_ fyl_dn 


ijtdt -|- adMf, 

\yt yt. 

J 2 

and 

have the 

following 

covariations 


M% = {yl - xl) dt, d{M\ My)t = -x^dt, d{M\ My)t = {yj - y^) dt, 
and can thus be represented by a real Brownian motion 

dM^ = ijtdBt, dMf = -xtdBf 

Recall moreover that the process is parametrized by the arc length so that we have almost 
surely, for all time t > 0: 


(A.2) 


+ y'i 

Vt 


= 1. 


Note that since the process can not explode in hnite time, we have yt G (0, +oo) for all 
f > 0 almost surely. As in the case where d > 3, the study of the process is then made easier 
thanks to the presence of two sub diffusions: 

Lemma A.1.1. The kinetic Brownian motion {xt,yt,Xt,ilt) admits the following 

subdiffusions of dimension one, two and three respectively: 


{ut)t>o ■— ( — 

-ht/ t>o 


(,yt,ut) 


t>o, 


{ytj Tt, yt)t>o- 


Proof. Starting from the system (A.l), a direct application of Ito’s formula gives 


dyt Uj^y^dt, 


dut = —(1 — u‘f)dt ——Utdt + cta/I — u^dBt, 
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hence the result for the hrst two sub diffusions. Moreover, from Equation (A.l), the evolution 
of {yt,Xt,yt)t>o does not depend on xt- □ 


Remark A.1.2. Note that, by symmetry if {yt,Xt,yt) is a version of the dimensional sub¬ 
diffusion starting from {yo,xo = 0,yo = ±1), then [yt, —Xt,yt) is also a version of the process. 


A.2. Long time asymptotics. From the ergodicity of the subdiffusion {ut)t>o, one can then 
deduce the following almost sure long time asymptotic behavior for the whole process: 


Proposition A.2.1. Let {xt,yt,Xt,ilt) be the kinetic Brownian motion in starting from 

{xo,yo,Xo,yo). Then the process {ut)t>o is ergodic in [—1,1] with invariant probability measure 

2 g-2x/(T2 

= zTrTli' 

where Z is a normalizing constant. Moreover, as t goes to infinity, we have 

lim -log ^ = [ xfi{dx) G (—oo,0), 

t^+oo t \yoJ J-i 

in particular y^ goes almost surely to zero exponentially fast and the process Xt converges almost 
surely to a random variable Xoo £ 1^- 



Figure 7. Simulation of the process and the corresponding geodesic for a = 0.01,0.1,1, 2. 


Proof. One easily checks that /i is an invariant probability measure for {ut)t>o which is thus 
ergodic. The ergodic theorem then ensures that almost surely 




Ugds / xpL{dx) G (—00,0). 


'0 J-i 

In particular, almost surely, yt goes to zero exponentially fast as t goes to inhnity. Moreover, 


from the relation (A.2), we have xf < yf from which we deduce that 

Xt = Xq = / Xgds 


is also almost surely convergent, hence the result. 


□ 
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A.3. Poisson boundary. We conclude the study of the kinetic Brownian motion in the 
hyperbolic plane by expliciting its Poisson boundary. We show that the invariant sigma held 
of the full process coincides almost surely the sigma held generated by the point (xoo, 0) G 
the boundary at inhnity. 

Proposition A.3.1. Let {xt,yt,Xt,iit) be the kinetic Brownian motion in starting from 

{xo,yo,xo,yo)- Then the Poisson boundary of the process coincides almost surely with a{xoo)- 


Proof. The proof is again an application of the devissage method. Since {ut)t>o is ergodic 
in [—1,1], its invariant sigma held is clearly trivial. Otherwise, the invariant sigma held of 
the process {yt,yt)t>o coincides almost surely with the one of the process {—\og{yt),Ut)t>o, 
whose behavior is stricly similar to the one of the process (r^, rt)t>Q in the proof of Proposition 
3.6.1, namely — log(|/t) goes almost surely to inhnity with t while its derivative Ut is ergodic. 


The exact same proof then shows that the invariant sigma held of {-~\og{yt),Ut)t>o and thus 
{yt,yt)t>o is trivial. Equivalently, if {yl,yl)t>o and {y^,y^)t>o are two versions of the process 
starting from (i/o,|/o) and (i/o)2/o) respectively, there exists two random times Ti and T 2 that 
are hnite almost surely and such that = (i/tj+o 2/T2+t)i>o- Without loss of 

generality, we can suppose that Ti> S where S is the hrst hitting point of ±1 by uj = yj/yj, 
which is also hnite almost surely since ui ergodic in [—1,1]. Now, the relation (A.2) implies 
that a shift coupling between {yj, yj) and {y'f, yl) yields automatically a shift coupling between 
the corresponding versions of the process {yt,Xt,yt,)t>o- Indeed, Equation (A.2) then ensures 
that \xff\^ = \xfff almost surely. If xf^ = xf^ , there is nothing to prove, if x^~= we can 

use the symmetry argument of Remark A.1.2 and consider the process (—starting at 
—Xg = 0. In other words, the invariant sigma held of the subdihusion {yt,yt,Xt)t>o is almost 
surely trivial. To apply the de vissage scheme, we are left to check the equivariance and 
regularity properties of Theorem 


3.2.8 


We have Xt = xq + fg Xgds so that the inhnitesimal 
generator of the full dihusion {xt,yt,Xt,yt) is clearly equivariant under the action of M by 
translation on the x component. Moreover, the generator is hypoelliptic so that harmonic 
function are continuous. By Theorem 3.2.8| , we can conclude that the invariant sigma held 
of the kinetic Brownian motion {xt,yt,Xt,yt) coincides almost surely with the sigma held 
generated by the limit Xoo- □ 
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